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Light scalar fields such as axions and string moduli can play an important role in early-universe 
cosmology. However, many factors can significantly impact their late-time cosmological abundances. 
For example, in cases where the potentials for these fields are generated dynamically — such as 
during cosmological mass-generating phase transitions — the duration of the time interval required 
for these potentials to fully develop can have significant repercussions. Likewise, in scenarios with 
multiple scalars, mixing amongst the fields can also give rise to an effective timescale that modifies 
the resulting late-time abundances. Previous studies have focused on the effects of either the first 
or the second timescale in isolation. In this paper, by contrast, we examine the new features that 
arise from the interplay between these two timescales when both mixing and time-dependent phase 
transitions are introduced together. First, we find that the effects of these timescales can conspire 
to alter not only the total late-time abundance of the system — often by many orders of magnitude 
— but also its distribution across the different fields. Second, we find that these effects can produce 
large parametric resonances which render the energy densities of the fields highly sensitive to the 
degree of mixing as well as the duration of the time interval over which the phase transition unfolds. 
Finally, we find that these effects can even give rise to a “re-overdamping” phenomenon which causes 
the total energy density of the system to behave in novel ways that differ from those exhibited by 
pure dark matter or vacuum energy. All of these features therefore give rise to new possibilities for 
early-universe phenomenology and cosmological evolution. They also highlight the importance of 
taking into account the time dependence associated with phase transitions in cosmological settings. 


I. INTRODUCTION 


Light scalar fields are a common feature in many 
cosmological scenarios and thus frequently play an im¬ 
portant role in early-universe cosmology. Such scalars 
include, for example, the QCD axion [1-4] and other 
axion-like particles [5], supersymmetric partner particles 
such as sneutrinos and staus [6-9], string moduli such 
as the dilaton and other geometric moduli [10-12], Q- 
balls [13, 14], quintessence fields [15] and chameleons [16- 
18], familons [19] and Majorons [20, 21], certain degrees 
of freedom within inert Higgs-doublet models [22], ad¬ 
ditional scalars present in little-Higgs theories [23], and 
others [24-27]. Some of these fields may even serve as 
candidates for dark matter or as significant contributors 
to dark energy. Moreover, the number of such scalars can 
be relatively large; indeed, the recent Dynamical Dark 
Matter framework [28, 29] posits a dark sector containing 
an entire ensemble of such fields which together conspire 
to produce a number of phenomenological, astrophysi- 
cal, and cosmological effects which differ markedly from 
those arising from more traditional dark sectors. Like¬ 
wise, string-inspired models typically involve significant 
numbers of moduli — and frequently a large number of 
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light axions and axion-like fields as well [30-33]. 

When a scalar field is light, the reason is usually that 
a mass term for this field is forbidden by a symme¬ 
try of the full, high-temperature theory. Such a mass 
term must therefore be generated at a lower tempera¬ 
ture scale through some dynamics which violates these 
symmetries — dynamics typically associated with a cos¬ 
mological phase transition. For example, masses for ax¬ 
ions or axion-like particles can be generated by non- 
perturbative instanton effects which become significant 
at temperatures near the confining temperature of the 
corresponding gauge group. Indeed, dynamical mass gen¬ 
eration through cosmological phase transitions provides 
a natural mechanism for engineering mass-scale hierar¬ 
chies. 

However, dynamical mass generation can also have po¬ 
tentially significant ramifications for the cosmology of the 
particles involved. First, cosmological phase transitions 
are never instantaneous. In the case of a second-order 
phase transition, the transition between the high- and 
low-temperature phases is smooth but nevertheless time- 
dependent. Even in the case of a first-order transition, a 
rapid epoch of bubble nucleation and expansion interpo¬ 
lates between the two phases. As a result, scalar masses 
which are dynamically generated during either kind of 
phase transition do not “turn on” instantaneously every¬ 
where in the universe; rather, these masses depend non- 
trivially on time and evolve continuously over the course 
of the phase transition towards whatever asymptotic val¬ 
ues they will have at late times. It turns out that this 
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time-dependence — and the existence of an associated 
timescale or “width” over which the phase transition un¬ 
folds — can have a significant impact on the cosmology of 
light scalars. For example, in the case of an axion receiv¬ 
ing its mass via non-perturbative instanton effects, the 
non-trivial consequences of this time-dependence on the 
resulting axion cosmological abundance were assessed in 
Ref. [34]. 

Another important consequence of dynamical mass 
generation is that it can give rise to mixing amongst 
scalars sharing the same quantum numbers. Indeed, if 
multiple such scalars couple to the same source of dynam¬ 
ical mass generation, there is no guarantee that the mass 
eigenstates of the theory will be the same before and af¬ 
ter the corresponding phase transition. The phenomeno¬ 
logical consequences of such mixing were investigated in 
Refs. [29, 35, 36] for an ensemble of light axion-like fields. 
Mixing can also be induced amongst light scalar fields at 
low energies through a variety of other dynamical pro¬ 
cesses, including the integration out of massive helds to 
which our light scalars are coupled. For example, the 
phenomenological implications of such mixings amongst 
axion-like fields have been studied in Ref. [37]. Regard¬ 
less of the source, however, we find that dynamically- 
induced mixing generally results in a non-trivial modi¬ 
fication of the coupling structure of the fields involved 
and can thus have a significant impact on the resulting 
energy densities of the fields at late times. 

The effects of mixing and of a non-trivial time- 
dependence for scalar mass generation have historically 
been studied separately. In this paper, by contrast, we 
shall examine the effects that arise when the mass matrix 
for a set of multiple light scalar fields includes both mixing 
and a non-trivial time-dependence. As we shall demon¬ 
strate, the interplay between the different timescales as¬ 
sociated with these effects can, in the presence of a time- 
dependent background cosmology, give rise to a number 
of qualitatively new features which are not seen when 
these effects are each considered in isolation. 

These new features generally concern the late-time cos¬ 
mological abundances (energy densities) associated with 
the scalar fields in our coupled multi-scalar system. For 
example, we shall demonstrate that the total late-time 
energy density of our system can be altered — often by 
many orders of magnitude — compared with traditional 
expectations. This includes situations in which the late¬ 
time energy density is enhanced to greater values as well 
as situations in which it is strongly suppressed. Second, 
we shall demonstrate that not only can the total energy 
density be altered, but so can the distribution of this to¬ 
tal energy density across the different fields in our system 
— often in dramatic ways. Indeed, these effects can even 
completely change which held carries the largest abun¬ 
dance. Third, we shall demonstrate that under certain 
circumstances, the combined effects of our two timescales 
can give rise to large parametric resonances which ren¬ 
der the energy densities of the helds highly sensitive to 
the degree of mixing as well as the duration of the time 


interval over which the phase transition unfolds. As a 
result of this heightened sensitivity, any effects which 
modify the mixing slightly or which change the rate at 
which our mass-generating phase transition unfolds can 
have a huge inhuence on the corresponding late-time 
cosmological abundances. Finally, we shall demonstrate 
that mixing in conjunction with a time-dependent phase 
transition can even give rise to a “re-overdamping” phe¬ 
nomenon which causes the field values and total energy 
density of the system to behave in novel ways which differ 
from those normally associated with pure dark matter or 
vacuum energy. 

Taken together, we see that all of these features give 
rise to new possibilities for the phenomenology of the 
early universe and cosmological evolution. In particu¬ 
lar, many new possibilities for model-building emerge. 
However, our results can also serve as warning: in the 
presence of multiple mixed scalar helds, it is essential to 
treat mass-generating phase transitions rigorously, with 
the proper time-dependence included. Indeed, approxi¬ 
mations in which such phase transitions are treated as in¬ 
stantaneous can produce late-time energy densities which 
differ from their true values by many orders of magni¬ 
tude. This warning is particularly relevant in the case 
of axions, where the mass-generating phase transition 
is nothing but the instanton-induced QCD phase tran¬ 
sition. This phase transition unfolds over a calculable 
non-zero timescale, and we shall demonstrate explicitly 
that ignoring this timescale (by setting it to be either 
zero or effectively inhnite) can lead to highly inaccurate 
late-time axion abundances. 

This paper is devoted to a general study of the new fea¬ 
tures discussed above. Because these features are com¬ 
mon to many systems of scalars which exhibit mixing 
in conjunction with a time-dependent mass-generating 
phase transition, in this paper we shall work within the 
context of a general model of scalars <j)i whose identi¬ 
ties remain unspecified. We shall likewise make no as¬ 
sumptions about the nature of the mass-generating phase 
transition except that it unfolds over a certain timescale. 
Moreover, quite remarkably, we shall find that a simple 
toy model consisting of only two scalar fields is sufficient 
to illustrate all of the features outlined above. This paper 
will therefore focus on an analysis of this two-component 
toy model, and we shall defer a study of more complex 
scenarios to future work [38] . 

This paper is organized as follows. In Sect. II, we intro¬ 
duce the two-component toy model which forms the basis 
for all subsequent discussions in this paper. As outlined 
above, this toy model contains both a mixing between our 
two helds as well as a time-dependent mass-generating 
phase transition. Along the way we also introduce all 
needed dehnitions, conventions, and notation. We also 
discuss some of the basic properties of this model. Then, 
in Sect. Ill, we study the behavior of the total late-time 
energy density of our two-component system, while in 
Sect. IV we study how this total late-time energy density 
is distributed between the two helds of our model. In 
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Sect. V we discuss the parametric resonance which ap¬ 
pears in the energy density, while in Sect. VI we turn our 
attention to the “re-overdamping” phenomenon. 

Although our analysis up to this point is completely 
general, our results have many immediate phenomeno¬ 
logical applications. In Sect. VII, we therefore focus on 
a special case by considering the effects on the standard 
QCD axion that emerge when a second axion is incorpo¬ 
rated into the theory. In Sect. VIII we then sketch some 
of the additional features that emerge when more than 
two fields are considered. These features will ultimately 
be explored more fully in Refs. [38]. Finally, in Sect. IX, 
we conclude with a discussion of our main results and 
possible avenues for extension and generalization. 

This paper also contains two Appendices. Although 
most of the results of this paper are obtained through 
numerical analysis. Appendix A provides exact analytical 
results in certain tractable special cases. Appendix B 
then discusses an alternative approach towards analyzing 
some of the features presented in this paper. 


II. A TOY MODEL 


In this section we delineate the toy model which shall 
be the basis of our analysis in this paper. This section 
will also serve to establish our notation and conventions. 


A. Motivation: A one-component warm-up 


As a prelude to the presentation of our model, we be¬ 
gin with a short discussion of the one-component case. 
Along the way we will also review some basic facts about 
scalar fields in an FRW cosmology and provide motiva¬ 
tion for the particular construction of our eventual multi- 
component toy model. 

Towards this end, let us consider a single scalar field (j) 
of mass m evolving in a flat FRW universe. If all spatial 
variations in (f> are assumed to be negligible, such a field 
evolves according to the standard equation of motion 

^ + 3i?(t)|^+m20 = 0 (2.1) 

where the Hubble parameter H{t) scales as 3H{t) njt 
with K = 2 (respectively k = 3/2) during a matter- 
(radiation-) dominated epoch. As a result, the field 
behaves as a damped oscillator with a time-dependent 
damping ratio Cit) = iH{t)/2m. Since H{t) falls with 
time, a field which is initially overdamped will inevitably 
become underdamped and experience decaying oscilla¬ 
tions. 

For simplicity, the energy-momentum tensor for such a 
field can be modeled as a perfect fluid with energy density 


p and pressure p given by 


(dip/dt)^ + 
{d(j)/dt)^ — 


( 2 . 2 ) 


If (j) is initially in an overdamped phase with d(j)/dt = 0, 
then p = —p: the energy density p associated with such 
a field in this phase behaves as vacuum energy. By con¬ 
trast, after (j) transitions to an underdamped phase, the 
field eventually experiences oscillations which are approx¬ 
imately virialized: in this phase p = 0 and the corre¬ 
sponding energy density p can be associated with massive 
matter. Of course, there is also an intermediate time in¬ 
terval during the transition from the overdamped to un¬ 
derdamped phase within which the behavior of the field 
exhibits transient features that eventually dissipate. 

The equation of motion (2.1) can be simplified by defin¬ 
ing a dimensionless time variable t = mt. The solutions 
can then be expressed analytically in terms of Bessel 
functions of the first and second kind. The exact evo¬ 
lution of the corresponding energy density p is shown 
in Fig. 1, where tq denotes the critical damping time 
at which C = 1- The different phases of the system are 
clearly evident in Fig. 1. For early times r ^ r^, the 
system is overdamped and the energy density p is essen¬ 
tially constant: p ~ po- For late times r ^ r^, by con¬ 
trast, the system is underdamped and the energy density 
p scales as p ^ Finally, during the intermediate 

times T ^ T^, the energy density p exhibits a correspond¬ 
ing transition between the two limiting behaviors above, 
punctuated by small transient oscillations. Although the 
specific dynamics of the field during this transitional pe¬ 
riod does not affect the late-time asymptotic scaling of 
p with time, these transients are nevertheless important 
in determining the overall scale of the late-time energy 
density. 

In this example, we have taken the mass to to be a 
constant, non-zero for all time. However, in many cosmo¬ 
logical situations, masses are generated by phase transi¬ 
tions. In such cases, the masses of such fields can be time- 
dependent. A well-motivated example of this is case of 
the QCD axion: the axion potential is flat (i.e., to = 0) at 
temperatures T ^ Aqcd, but this flat potential is mod¬ 
ified at lower temperatures by instanton effects which 
generate an effective mass. 

Time-dependent masses to(t) can signihcantly modify 
the evolution of such scalar helds. In this paper, we are 
primarily interested in cases in which at least a portion 
of the contributions to the scalar masses are generated as 
the result of phase transitions. We shall therefore gener¬ 
ally consider cases in which to(t) has one (smaller) value 
at early times, a second (larger) value at asymptotically 
late times, and a smooth time-dependent transition be¬ 
tween the two. Indeed, we shall let tq represent the 
“central” time at which this mass-generating phase tran¬ 
sition occurs, and imagine that this transition unfolds 
over a time interval of duration or width Aq. We then 
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t/tq 

FIG. 1: The energy density p, normalized to its initial value 
po, for a solution to Eq. (2.1) in a radiation-dominated (RD) 
universe, i.e., n = 3/2 (red), or a matter-dominated (MD) 
universe, i.e., k = 2 (blue). In the overdamped phase, the en¬ 
ergy density remains nearly constant at p = —p > 0, behav¬ 
ing as vacuum energy. By contrast, once the critical damp¬ 
ing C = 1 threshold is crossed at T( (a critical time which is 
slightly different for RD and MD universes), the energy den¬ 
sity begins to dissipate, asymptotically exhibiting the simple 
power-law scaling behavior p ~ expected for matter. 

obtain the situation illustrated schematically in Fig. 2. 

The introduction of a non-zero time interval Aq dur¬ 
ing which our phase transition unfolds can have a signif¬ 
icant effect on the time-evolution of (j). When Ac = 0, 
there are essentially two regimes which one can consider, 
each with its own distinctive phenomenology: < tq 

and tq > tq. However, a non-zero width opens up 
a third possibility, with ~ tq. Indeed, more gener¬ 
ally, the introduction of a non-zero Ac introduces a new 
timescale into the problem, and we shall see that this can 
have dramatic effects. We shall explore the phenomenol¬ 
ogy of this region as part of our larger general study. We 
emphasize, however, that phenomenon of having ~ tq 
is completely different from having tg during the “damp¬ 
ing transition” region of Fig. 1. Indeed, while Aq may be 
considered to represent a width around tq, the “damp¬ 
ing region” of Fig. 1 can instead be considered a width 
around r^, a width which exists — as shown in Fig. 1 — 
even if the mass is constant. 


B. Defining our toy model 

The introduction of non-zero time interval Aq for our 
phase transition is just one feature we wish to study in 
this paper. The other concerns the possibility of mix¬ 
ing between different held components tpi, i = 1,..., N. 



FIG. 2: Two time-dependent mass scales: a growing mass 
function m(r) and a falling Hubble friction scale 3H{t)/2. 
Their intersection determines the transition time p;. We as¬ 
sume that the growth of m(T) takes place during an interval 
of approximate duration Aq centered at tq. As we shall see, 
the introduction of a non-zero Aq can have many dramatic 
effects on the late-time energy densities of those fields which 
(either directly or indirectly) are sensitive to this change in 
mass. For example, one new possibility is that the transition 
between overdamped and underdamped phases can occur dur¬ 
ing the mass-generating phase transition. 


The phenomenological effects stemming from each can 
certainly be studied in isolation. However, as we shall 
hnd, these two features can conspire to produce a num¬ 
ber of remarkable effects that transcend what is possible 
with either alone. For this reason, we shall study the 
implications of both features together. The effects of ei¬ 
ther feature in isolation can then be extracted through 
the limiting cases in which the effects due to the other 
feature are gradually turned off. 

In this paper, we shall explore these effects within the 
context of a simple toy model involving only two compo¬ 
nents, (jjo and (fi. As we shall see, our toy model is simple 
enough to be tractable, yet rich enough to incorporate all 
of the phenomena of interest. 

Our model consists of two real scalar fields, 4>q and (fi. 
If we again assume that the spatial variations in these 
fields are negligible, their equations of motion in a flat 
FRW universe take the form 

3 

where \s the corresponding squared-mass matrix. At 
times t to, long before mass generation occurs, we 
shall take A4^ constant and diagonal. In fact, for simplic¬ 
ity, we shall further assume that (po is massless at such 
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early times. Thus, at early times, we shall assume 


- 

t<s:tG 


0 0 

0 M2 


(2.4) 


where M ^ 0 is a general unfixed mass parameter. By 
contrast, long after the phase transition has occurred, we 
assume that new components mT will have been gener¬ 
ated in the squared-mass matrix: 


- 

f^ta 


0 

0 

+ 


0 

m2 

Wqi ™11_ 


(2.5) 


In order to connect these two asymptotic extremes, we 
shall let 'rnfj{t) denote time-dependent elements of the 
mass matrix which interpolate between zero at early 
times and at late times. We can then write 


M^it) 


0 0 

0 M2 


mooit) ™oiW 

mhit) TO?i(t) 


( 2 . 6 ) 


Having specified the mass matrix, we can now intro¬ 
duce a dimensionless time variable t = Mt, as in the 
single-component case. Our equations of motion then 
take the form 

0,+3M(r)(^,+^M(2^-(r)0, = 0 (2.7) 

3 

where the dots indicate d/dr and where our mass matrix 
is also dimensionless and takes the form 


M^{t) 


WooM TOqiM 
TOgi(r) l-hmfi(r) 


( 2 . 8 ) 


with each dimensionless rriij now understood to be a frac¬ 
tion of M. We shall adopt these conventions throughout 
the rest of this paper. 

In general, there is no reason to expect that the mass 
generation occurs in the same basis as in Eq. (2.4). We 
shall therefore allow for the possibility that 0 — 

be., the possibility that our phase transition induces a 
mixing between our primordial field components (^o and 
(fi- However, we shall nevertheless make the simplify¬ 
ing assumption that the time-dependence of each com¬ 
ponent is identical, allowing us to focus on those effects 
that come from a single common timescale for mass gen¬ 
eration. Indeed, since the mass matrix is nothing but the 
curvature matrix associated with the potential V[4>q^4>i) 
induced by the phase transition, this assumption is tan¬ 
tamount to assuming a single time-dependence for the 
potential as a whole. As a result, we can write each of 
the individual mass components in a factorized form 


mijir) = niij ■ h{T) , (2.9) 

where ^(t) is a smooth function of time which describes 
the time-development of the phase transition. 


Our next step in specifying our toy model is to choose 
a suitable function /i(t). As indicated above, we require 
that /i(r) —0 as r —0 and /i(r) —1 as r —>■ oo. How¬ 
ever, because one of our main interests in this paper con¬ 
cerns the timescale associated with the mass-generating 
phase transition, we would also like /i(r) to incorporate a 
dimensionless parameter a which controls how abruptly 
the phase transition occurs. The limit cr —> 0 might then 
correspond to a phase transition which is effectively in¬ 
stantaneous, while non-zero values of cr correspond to 
phase transitions which occur increasingly slowly. It is 
also desirable, regardless of the width of the transition, 
that the midpoint at tq be a fixed point of reference. 
We therefore include in our construction the requirement 
that h{Tc) = 1/2 for all a. Indeed, this may be taken as 
a definition of tg- 

Beyond these constraints, the choice of /i(r) is com¬ 
pletely arbitrary, and many functions may be chosen. For 
concreteness, however, we shall take 


h{T) 


1 

2 


11 -|- erf 


- log 
cr 



where the error function erf(z) is given by 
erf( 2 ) = [ e~^ dx . 

Jo 


( 2 . 10 ) 


( 2 . 11 ) 


As illustrated in Fig. 3, this function satisfies all of our 
requirements. Of course, many other choices for h(r) are 
possible. However, none of the qualitative results of this 
paper will ultimately depend on the specific choice for 
/i(t). Thus, any smooth, monotonic function /i(t) satis¬ 
fying the above constraints will lead to similar results. 

Corresponding to each non-zero value of the parameter 
a there exists a non-zero timescale Aq over which the 
phase transition occurs. In general, we may define Aq 
in terms of the slope of the hir) function at its midpoint 

tg- 


h{TG + St) « o + for (5r <C Ag . (2.12) 

2 Ag 

Adopting this definition for Ag, we then find for our /i(r) 
function that 


Ag = VraTG ■ (2.13) 

We thus see that Ag —0 as ct —0, as expected. Indeed, 
this limit corresponds to the case of an instantaneous 
phase transition, with the hir) taking the form of the 
Heaviside step function 0(r — tg). 

In this paper, we shall consider the effects that arise 
when cr is non-zero. There is, however, a critical value 
a* above which the behavior of /i(t) changes in an im¬ 
portant way. To see this, let us first consider the limit 
in which ct —>■ oo. In this limit, we find that Ii(t) is es¬ 
sentially constant and rG-independent: /i(t) « |0(t) for 
all r. This limit may therefore be interpreted as one in 
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FIG. 3: The general form of the h-function in Eq. (2.10), 
plotted for several choices of the parameter a. In all cases we 
see that /i(t) rises smoothly from h — 0 to h = 1 and crosses 
h — lj2 at T = TG- Large values of a correspond to more 
gradual transitions, while for cr — >■ 0 we find that h approaches 
the Heaviside step function 0(t — tq)- 


which our original phase transition at tq effectively dis¬ 
appears and is replaced by a new, infinitely sharp “phase 
transition” at t = 0. Note that this latter transition is 
nothing but an artifact of our boundary condition that 
/i = 0 at T = 0. As a result, dialing a from 0 to oo has 
the effect of slowing our original phase transition near tq 
while simultaneously building up an artificial phase tran¬ 
sition near r = 0. Indeed, if a grows too large, the most 
rapid changes in the mass parameters of our system will 
no longer be associated with our original phase transition 
near r = tg, but with the artificial one near r = 0. 

In this paper, we wish to maintain the notion that 
increasing the value of cr corresponds to slowing the time- 
development of our mass matrix. Even more importantly, 
we also wish to ensure that the time period exhibiting the 
most rapid time-development of our mass matrix is still 
associated with our original phase transition near tq and 
not the artifact near r = 0. Therefore, in this paper, we 
shall always restrict our attention to values 0 < tr < ct* 
where a* is that value of a for which maXr{dh/dT) is 
minimized. For the hir) function given in Eq. (2.10), we 
find that cr* = for all rg. As a result, in what follows 
we shall only consider values of Aq in the range 

0 < Ag < V^tg ■ (2.14) 

In order to gain insight into the values of a that we 
might expect for a well-motivated phase transition, let 
us consider the case of the instanton-induced phase tran¬ 
sition during which a mass is generated for the QCD 
axion. Detailed lattice studies [39, 40] have yielded ap¬ 


FIG. 4: The function /iqcd describing the instanton-induced 
time-dependent mass of the axion (solid blue line), plot¬ 
ted as a function of t/ta- This curve is generated for 
Aqcd = 200 MeV, with to defined as that location where 
the curve crosses h = 1/2. For comparison purposes, our h- 
function in Eq. (2.10) with a ~ 0.47 is superimposed (dashed 
red line). We see that our h-function matches /iqcd surpris¬ 
ingly well throughout the entire range of values of plotted, and 
fits particularly well near the “center” of the phase transition 
at t = to • 


proximate expressions for the mass of the QCD axion 
as a function of temperature T. However, the cosmologi¬ 
cal time-temperature relation appropriate for a radiation- 
dominated epoch is given by 



where Mp is the Planck mass and where g* is the 
temperature-dependent effective number of relativistic 
degrees of freedom. Use of this relation then allows us to 
determine the axion mass as a function of cosmological 
time t, with the results shown in Fig. 4. For comparison, 
our h-function with a « 0.47 is also shown in Fig. 4. It is 
clear that these two functions match quite well over the 
entire range of times shown. 

The only remaining ingredient to be specified as part of 
our toy model is a set of initial conditions to be imposed 
at some early time tq. For the differential equations in 
Eq. (2.7), this means specifying initial values for our two 
fields 00,1 and their first derivatives 0o,i- While many 
possibilities exist, one particularly natural choice is to 
consider a mere displacement for the 0o field: 


( 0o\ _ ( Ao\ ( 00 

1^01/ " I 0 j ’ i0i 


0 

0 


(2.16) 


r=To 


r=To 
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There are three special properties associated with this 
initial configuration which make this choice especially ap¬ 
pealing. First, as long as tq ^ tc, this configuration does 
not introduce any initial energy into the system. Thus, 
all of the energy that our system accrues will be solely 
that injected through the phase transition. Second, as 
long as Tq Tq, the specific choice of the initial time 
To will be irrelevant. In other words, given our other as¬ 
sumptions, the system is essentially time-independent for 
all T TG and thus completely insensitive to variations 
of Tq. Asa result, we shall never need to refer to tq again, 
understanding implicitly throughout this paper that tq 
has always been chosen to be sufficiently early that our 
results are independent of tq. Finally, the choice of Aq 
in Eq. (2.16) merely serves to set an overall mass scale 
for our system. However, since we shall always be con¬ 
sidering ratios of field values or energy densities in this 
paper, our results will also be insensitive to Aq. 

This completes the specification of our toy model. It 
contains only five free parameters: two (namely tq and 
Ag ^ a) which describe the time-development of the 
phase transition, and three (namely TOqq, and 
which describe the effects of the phase transition and 
the mixing it induces. Despite its simplicity, however, we 
shall find that this toy model is not only of sufficient gen¬ 
erality to accommodate a wide variety of physical systems 
but also of sufficient richness to give rise to a number of 
surprising phenomena. 


C. Preliminary analysis: Constraints on mixing 


However, for many purposes it will prove useful to intro¬ 
duce the variables 


= Mn+Moo = l + mlo + mj^ 

Am^ = MIi-MIq = l-mlQ + mh (2.19) 

in terms of which our squared-mass matrix in 

Eq. (2.8) takes the more symmetric form 




5 ("iLm - wgi 

Woi 5 {rrisum + 


( 2 . 20 ) 


Note while the individual mfj quantities each carry the 
same time-dependence [proportional to /i(t)^], the same 
is no longer true for and Arn^. In terms of these 

variables, the constraints in Eq. (2.17) then take the form 


[ • > 0 

< • [Am^l < _ (2.21) 

[• [moil < - (Am2)2 , 

whereupon we find from Eq. (2.18) that 

[wgilmax = - {Am?f . (2.22) 

Finally, perhaps the most useful way to parametrize 
the mixing in our toy model is in terms of a rotation 
angle 9 which relates the mass eigenstates (j)\o (j)\^ 

at any instant of time to the original mass eigenstates 
prior to the onset of the phase transition: 


In subsequent sections, we shall analyze the behavior 
of our fields and of their corresponding energy densities 
within different regions of our five-dimensional parameter 
space. However, even before proceeding, there are a num¬ 
ber of preliminary observations that hold more generally 
and which serve to significantly constrain the allowed pa¬ 
rameter space. 

These constraints all ultimately stem from the obser¬ 
vation that not every choice of the masses can be 
made independently. Requiring that the eigenvalues of 
Al^ be real is tantamount to demanding that be 
Hermitian. However we must further demand that these 
eigenvalues be non-negative, which requires that Al^ be 
positive-semidefinite. This requires that all be real 
and satisfy the three constraints 

[ • rnlo>0 

< • mfi > —1 (2-17) 

[• wgo(l + "iii) > ^01 ■ 

These constraints provide a set of intrinsic limits on 
the magnitude of the mixing that may occur within our 
toy model. Indeed, these constraints imply that 

|tooi| < [w^ilmax = \/"roo(l + "^ii) ' (2-18) 


_ f(^ose -sin6»\ f(j)o\ 

[sine cose) [cj),) ’ 

from which it follows that 

tan(26') = . (2.24) 

Am^ 

Indeed, while 0 generally populates the range 
—TT < 9 < n in Eq. (2.23), we see from Eq. (2.24) 
that it is sufficient to focus on —tt/ 2 < 9 < Trj2 for the 
purpose of calculating mixing angles, since these are 
invariant under the mapping ^ 

The same will also be true for calculating energy 
densities — our main interest in this paper — because 
these energy densities will depend only quadratically on 
the fields in Eq. (2.23). However, we further observe 
that in the ^Ai} basis introduced in Eq. (2.23), the 
equations of motion in Eq. (2.7) now take the form 

^Ao + 3i7^Ao+ (-^0 ~ '(’Ao 

= - 20,^Ai - (0 + 3770) 

4>\i + 377(^Ai+ (^^1 — ^^) (f’Xi 

= + 29^X0 + {d + 3770) c^Xo ■ ( 2 . 25 ) 








These equations of motion exhibit an invariance 
under the simultaneous correlated transformations 
{4>Xg, and 9 —9. Since this is also 

an invariance of our initial conditions in Eq. (2.16), our 
toy model will exhibit this invariance at all points during 
its time-evolution. Note that this invariance is also tan¬ 
tamount to the simultaneous correlated transformations 
(d’Oj'/'i) and rri^i —>■ —TOq]^. We can hence¬ 

forth truncate our attention to rri^i > 0 and 0 < 0 < 7r/2 
without loss of generality, and we shall do so through¬ 
out the rest of this paper. The region 0 < 0 < Tr/d then 
corresponds to > 0, while the region 7r/4 < 6* < 7r/2 
corresponds to < 0. 

It then follows from Eq. (2.24) that for any values of 
and Am^, the allowed mixing angles 9 are those 
for which 

|tan(20)| < |tan(20,,ax)| = -1 (2.26) 

where 0max is defined as the value of 9 which maximizes 
I tan(20)|. Thus 9 < 6>max for 6<max < 7r/4, while 9 > 0max 
for 0max > Tr/d. These angles are illustrated in Fig. 5. 
Note that for Am^ = (corresponding to Wqq = 0 

or mil = ”1) respectively), the corresponding values of 
9 are restricted to 0 or 7r/2. Indeed, any other values 
would lead to a squared-mass matrix with at least 
one negative eigenvalue, corresponding to a tachyonic 
mode. Finally, for Am^ = 0, we have 9 = tt/A for all 
mQi > 0, while the angle 9 is in principle undetermined 
for mQi = 0. 

We conclude this section with some further definitions 
that shall also prove useful in what follows. Quite often, 
we shall be evaluating a quantity X which is a function 
of the time r. For example, such quantities might include 
the mass parameters TOy(T), the mixing angle 9{t), the 
field values 4'xi'r), or the total energy density p ( t ). For 
any such quantity X{t), we shall define 

A = lim X{t) . (2.27) 

T—^OO 

In other words, X shall denote the asymptotic late-time 
value of A(r), where in practical terms the notion of a 
“late” time can be taken as referring to a time at which 
both fields (f>x have reached the asymptotic underdamped 
regime, with corresponding energy densities px exhibiting 
the virialized damped scaling behavior px ^ shown 
in Fig. I. Finally, we shall also find that a crucial mea¬ 
sure for the degree of mixing that might be present in 
a given system is not the absolute value of the mixing 
parameter m^i (or 0), but rather the value of this pa¬ 
rameter as a fraction of the total degree of mixing that 
would have been allowed for that system, given the con¬ 
straints discussed above. The same is often true for the 
splitting parameter Xrnf. Towards this end, we define 
the mixing saturation 

c = ^01 ^ |tan(26>)| ^ m^ 

[”^oi]max |tan( 26 >max)| y^mgo(I -|- mh) 

(2.28) 



FIG. 5: Allowed ranges for the mixing angle 9 (blue shaded 
region), plotted as a function of Am'^/rnsum- Note that 
9 is restricted to 0 or 7r/2 when Am^ = irUgum — i.e., 
when moo = 0 or mfi = — 1, respectively. By contrast, when 
Am^ = 0, the allowed range for 9 depends on the value of 
m§i: for moi A 0, the only allowed 6-value is 6 = 7r/4 (blue 
dot), while for moi = 0 we find that tan(26) is indeterminate 
and any value for 9 is allowed (red line), depending on how 
relevant limits are taken. 


as well as the analogous splitting saturation 

_ Am? Am? 1 — m^Q mh 

^ mhrn 1 + "*oo + "^11 

Note that in terms of these variables^ we have 
tan(20) = - \/l — 77^ . 

77 


( 2 . 29 ) 


( 2 . 30 ) 


We emphasize, however, that a maximally saturated mix¬ 
ing configuration with ^ = 1 does not necessarily imply 
that the mixing itself is maximal or even large on an 
absolute scale. For this reason, we shall usually specify 
the degree of mixing in a given configuration by quoting 
both ^ and ^max- Indeed, for many quantities, it will be ^ 


^ We caution the reader that unlike the other variables we 
have thus far introduced, neither ^ nor rj has a unique time- 
dependence. In other words, knowledge of the value of either ^ 
or 77 at a given time does not fix the value of ^ or 77 at other 
times without knowledge of the values of the more primordial 
variables from which ^ and 77 are derived. Thus, when we 
write a relation such as that in Eq. ( 2 . 30 ), we are illustrating 
the functional dependence of 0 on ^ and 77 and asserting that 
9 does not depend on the more primordial variables except in 
these combinations. However, the time-dependence of 9 cannot 
be determined from such an expression and must be calculated 
directly in terms of the primordial variables themselves. 
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rather than 9 or 0niax which characterizes the behavior of 
interest and which allows us to compare across systems 
with different values of 9 or 0niax- 

Note that for Am.^ = 0, ambiguities can arise when 
defining In particular, the middle expression in 
Eq. (2.28) is formally indeterminate in such situations. 
In this paper, for Am^ = 0, we shall therefore define ^ 
through the final expression in Eq. (2.28). Thus ^ can 
vary even though 9 may be fixed at 7r/4. Moreover, again 
following the final expression in Eq. (2.28), in this paper 
we shall define 

= i™ (^lAm2=o) (2-31) 

for any quantity X. Thus the ^ —>■ 0 limit will always be 
a smooth one with 9 = 7r/4, even for Am^ = 0. 


D. Temporal properties of the model: 

Basic features 

Given the definition of our toy model, it is now rel¬ 
atively straightforward to study the corresponding dy¬ 
namics. Indeed, this can be done numerically if not an¬ 
alytically, and certain features are entirely as expected 
and relatively easy to understand. However, many other 
features are surprising and will play a significant role in 
what follows. In the remainder of this section, therefore, 
we shall discuss general features of the time-evolution of 
this model. In particular, we shall focus on the time- 
dependence of the mass eigenvalues, mixing angles, and 
mass eigenstates — quantities upon which our future re¬ 
sults will rest. 

We begin by studying the mass eigenvalues in our 
model. At any moment in time, these masses are given 

by 

I =F yjtf + {I- 

= (1 Try sec 2d) . (2.32) 

The late-time values of these masses are shown in Fig. 6, 
while the evolution of these masses from early to late 
times is shown in Fig. 7. Numerous features are imme¬ 
diately apparent. First, as predicted from Eq. (2.32), we 
see that the eigenvalues at early and late times (r <C tq 
and T ^ tq, respectively) are independent of the sign of 
Xm? (or ry). [In this connection we recall from Fig. 5 that 
sec(2d) < 0 when rj < 0.] Thus, the early- and late-time 
values of Ag i are identical in the left and right panels of 
Fig. 7. However, we see that the time-evolution of the 
eigenvalues between these two endpoints is highly sensi¬ 
tive to the sign of Am^. For Am^ > 0, the eigenvalues 
evolve from initial to final values without any tendency 
towards level-crossing. For Am^ < 0, by contrast, the 
eigenvalues initially head towards each other as if to ex¬ 
perience a level-crossing. However, whether this level¬ 
crossing actually occurs depends on the value of the mix¬ 
ing. For ^ > 0, the non-zero mixing between mass eigen¬ 
states induces a level repulsion which ultimately prevents 


"^ 0,1 — 2 ^' 


1 ^2 

sum 





FIG. 6: The asymptotic (late-time) eigenvalues Xi plotted as 
functions of the late-time splitting and mixing parameters fj 
and G For all values of fj, we see that Ag —> 0 and Ai —>■ rfiaum 
as ^ ^ 1. 


a direct level-crossing. As a result, the eigenvalues veer 
away from each other, and ultimately assume the same 
late-time values that they had for Am^ > 0. 

The case with Afn^ < 0 and ^ = 0 is more subtle, and 
deserves special discussion. As ^ 0, we see from the left 

panel of Fig. 7 that our eigenvalues Ag and Ai actually 
meet near tq before bouncing off each other. [Note, in 
this connection, that Ag is always defined as the lighter 
eigenvalue, consistent with Eq. (2.32).] In this sense, no 
actual level-crossmyy occurs, even for ^ = 0; instead, each 
eigenvalue develops a “kink” — i.e., a discontinuous slope 
— at their meeting point near t = tq. 

Finally, for Am^ = 0, we see from Fig. 7 that there is 
no tendency towards level-crossing. However, for ^ = 0, 
the two eigenvalues approach each other asymptotically. 

We can also study the mixing angle 6* as a function of 
time. The result is shown in Fig. 8. As we see, for ^ > 0 
the mixing angle begins at zero, as expected, and rises 
to a non-zero late-time value 0; this corresponds to Am^ 
transitioning from its initial value Am^ = = I to 

its final value Afn^. Once again, however, the behavior of 
the angle 9 is highly sensitive to the sign of the late-time 
value Afn^. For Afn^ > 0, the angle 9 remains below 
7r/4 (i.e., within the lower right “lobe” of Fig. 5). For 
Am^ < 0, by contrast, the angle 9 eventually grows above 
7r/4, moving from the lower right lobe of Fig. 5 to the 
upper left. The case with Am^ = 0 will be discussed 
below. 

It is natural to define the time rg at which 9 crosses its 
midpoint value 9/2; in general, this is the time for which 

h(Te) = , \ = ■ (2.33) 

V1 + Ai - Ag 
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FIG. 7: The time-evolution of the two mass eigenvalues Aq^i as the phase transition unfolds, plotted for different mixing 
saturations. Typical behaviors are shown for Am^ < 0 (left panel), Am^ = 0 (middle panel), and Am^ > 0 (right panel). In 
all cases, the eigenvalues begin at {0,1}, as expected. For Am^ > 0, the eigenvalues slowly transition to their late-time values 
without any tendency towards level-crossing. For Am^ < 0, by contrast, the eigenvalues meet and rebound off each other in the 
case of zero mixing (^ = 0) but such a meeting is thwarted by level repulsion for all non-zero mixing. In general, the strength of 
the level repulsion increases with degree of mixing Despite these different features at intermediate times, we observe that the 
eigenvalues ultimately arrive at the same late-time values regardless of whether Am^ is positive or negative, in accordance with 
Eq. (2.32). Finally, for Am^ = 0, we see that there is no tendency towards level-repulsion, but for ^ = 0 the two eigenvalues 
approach each other asympotically. 


Likewise, for Am^ < 0, we find that the transition be¬ 
tween lobes at 9 = tt/A occurs at the time t* for which 

Hr*) = ■ ; . (2.34) 

V 1 — ArH 

In general, this “lobe-crossing” time r* is distinct from 
tq and Tg. However, t* and tq coincide when 9 = 7r/2. 
As we see from Fig. 8, this happens only when ArH < 0 
and ^ = 0. 

Once again, this case with Am^ < 0 and ^ = 0 
deserves special mention. At early times we have 
Am? = = 1, and our mixing angle 9 begins at zero. 

This angle then remains at zero all the way until the 
eigenvalue meeting time shown in the left panel of Fig. 7; 
note that this eigenvalue meeting time is indeed nothing 
but T* = Tg. At this time, 9 changes instantaneously to 
9 — 7r/2, consistent with the fact that after r* we have 
Am? < 0. Thus we see that 9 behaves as a step func¬ 
tion in this limit, with an effectively instantaneous field 
rotation. Note that this “instantaneous” behavior in the 
f —>■ 0 limit applies only for An? < 0. 

The situation with An? = 0 also deserves special at¬ 
tention. As usual, we have 0 = 0 at early times. For all 
^ 7 ^ 0, this mixing angle 9 then transitions to d = 7r/4 at 
late times. Moreover, as discussed above, it follows for¬ 
mally from Eq. (2.31) that 9 = 7r/4 even for ^ = 0. How¬ 
ever, there is an important subtlety in the latter case. As 
^ 0 with Am^ = 0, it turns out that h(rg) = 1. This 

in turn implies that Tg —>■ oo! Thus, even though 9 = tt 14, 
in this case, we never actually reach the point at which 
the corresponding mixing occurs. Rather, our two origi¬ 
nal states 00 and 0i remain unmixed at all finite times. 


Indeed, in some sense, the transition to 9 = tt/A occurs 
only when the two corresponding eigenvalues Aq and Ai 
actually meet. Note, in this connection, that even though 
the field rotation does not occur at any finite time, the 
mass eigenvalues Aq i nevertheless experience their nor¬ 
mal evolution in the neighborhood of tg, a result which 
follows directly from the time-dependence of the origi¬ 
nal mass matrix and which is independent of mixing. In 
other words, in this special case, the two timescales tq 
and Tg are maximally separated. 

Our discussion thus far has focused on the eigen¬ 
value/mixing structure of the mass matrix as a function 
of time. However, the corresponding field values also gen¬ 
erally behave as one might expect, at least as far as their 
grossest features are concerned. For times t <C rg, our 
two fields 0Ao 1 are effectively uncoupled: 4>\g evolves in¬ 
dependently of 0A^, and indeed 0 Ai remains vanishing. 
For times t ^ Tg, by contrast, the phase transition gen¬ 
erates a non-zero, time-dependent mixing which couples 
the two fields together and thereby causes 0 Ai to accrue 
a non-zero value as well. Finally, for times t ^ Tg, our 
mixing angle 9 is non-zero but essentially constant. This 
means that a similar decoupling exists during this period 
as well, except that our decoupled fields are now those lin¬ 
ear combinations which are rotated relative our original 
fields by 0. In general, each of these linear combinations 
will have a non-zero field value, and will evolve inde¬ 
pendently according to whether it is individually over¬ 
damped or underdamped. 

The limiting case with Am^ < 0 and ^ = 0 is again 
worthy of special note. In this case, 0 Ao retains its origi¬ 
nal amplitude until the eigenvalue meeting time r* = Tg. 
The entire amplitude of 0 Ao then transfers instanta- 
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n 



t/tg 


FIG. 8: The mixing angle 6, plotted as a function of time for 
Am^ = ±2 and different values of For early times t <^tg, 
we generally have Am^ = = 1 and 6 — 0. The subse¬ 

quent time dependence then depends critically on the sign 
of the late-time value Am^. If Am^ < 0, then Am^ eventu¬ 
ally switches from positive to negative values; this occurs at 
r = r*, where r* is defined in Eq. (2.34). At this point, the 
mixing angle 6 passes from the lower “lobe” of Fig. 5 to the 
upper lobe through the only allowed transition point between 
the two lobes at Am^ = 0 and 9 = rr/l. For Am^ > 0, by 
contrast, the mixing angle stays entirely within the lower lobe 
of Fig. 5 and thus remains below 7r/4. Finally, for Am^ = 0, 
the angle 9 moves to the central point between the two lobes, 
consistent with our expectation that 0 = 7r/4. However, as 
discussed in the text, it takes infinite time to do this as ^ 0. 


neously to (j)\i i where it remains for all later times. (This 
follows from the observation that in this limit we have 
an instantaneous rotation of our fields from 0 = 0 to 
0 = 7r/2.) This rapid amplitude transfer is consistent 
with the instantaneous infinite value of 9 at t* = rg. As 
a result of this amplitude transfer, quantities which de¬ 
pend on the field amplitudes (such as the associated en¬ 
ergy densities) transfer instantaneously from (j)\g to 
at this time. In other words, they “ride” smoothly across 
the eigenvalue collision in the left panel of Fig. 7, transi¬ 
tioning from (f)\^ to (^Ai ■ Of course, this is completely as 
expected, reflecting nothing more than the fact that the 
energy density began in our original field (j)^ and remains 
there when ^ = 0. All that has suddenly changed at the 
eigenvalue collision time is the identification between this 
field and our mass-eigenstate fields 4>\^^ ^. 

One could, in principle, continue along this line of 
inquiry. For example, one could proceed to study the 
phase-space trajectories of our mass-eigenstate fields ^Ao 
and (j)Xi 1 and map out how these trajectories depend on 
the different defining parameters of our model. However, 
an exhaustive study of this toy model is not our purpose 


in this paper. Rather, as stated in the Introduction, in 
this paper our interest in this toy model stems from the 
fact that — despite its simplicity — it gives rise to cer¬ 
tain features which have the potential to transcend our 
typical expectations when interpreted in a cosmological 
setting. It is therefore to these new features that we now 
turn. 

In keeping with the above observations, in the rest 
of this paper we shall mostly concentrate on those re¬ 
gions of parameter space in which our fields are already 
underdamped — or are in the process of becoming un¬ 
derdamped — during the mass-generating phase tran¬ 
sition. This is important, since the time-dependent ef¬ 
fects of our mass-generating phase transition tend to be 
washed out if our fields remain overdamped while it oc¬ 
curs. Thus, in this way, we shall be focusing on precisely 
those parameter-space regions of interest: those which 
are likely to exhibit a non-trivial interplay between the 
width of the mass-generating phase transition, the tran¬ 
sition between overdamped and underdamped regimes, 
and the mixing between all of the fields experiencing 
these effects. 


III. TOTAL LATE-TIME ENERGY DENSITY 


The total energy density p is the quantity of central 
interest in this paper. The role it plays as a cosmologi¬ 
cal observable gives it direct importance in any analysis 
of our toy model — particularly at late times, after our 
phase transition has been completed. In general, we are 
interested not only in the total late-time energy density 
p, but also in its distribution between the two individual 
components px- Moreover, in each case, we are particu¬ 
larly interested in knowing the extent to which our mass¬ 
generating phase transition might leave imprints on these 
late-time energy densities. These are the issues that we 
shall study in this section. 

Calculating the energy density of the system proceeds 
directly from the equations of motion for our two fields 
and their derivatives. In the original {(j)o,(pi} basis, the 
total energy density is given by 




(3.1) 


By contrast, in the mass-eigenstate basis {4>Xoi4>\i} 
troduced in Eq. (2.23), the total energy density is given 
by 



^ [4,1 + ( a > + 




XX' 


(3.2) 

where exx' is the Levi-Civita symbol with caqAi = +1- 
The first thing we observe is that in neither case can 
we express our total energy density as a sum of individ¬ 
ual contributions. Indeed, we cannot write p in the form 
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p = Pi ov p = Px- The reason for this, as most 
evident from Eq. (3.2), is ultimately that our fields ex¬ 
perience a time-dependent mixing, with 0^0. Without 
mixing — or even with only a constant mixing — such a 
decomposition could have been done and individual con¬ 
tributions identified. Thus, in this paper we shall never 
refer to the individual contributions to the total energy 
density except at late times when the mixing has essen¬ 
tially stabilized and 0 = 0. Indeed, at late times, we can 
then identify 

Px = 2 

where the overbar indicates late-time values in accor¬ 
dance with Eq. (2.27). 

As already discussed at the end of Sect. II, the effects 
of the mass-generating phase transition are washed out 
if our fields are still overdamped when it occurs. Like¬ 
wise, given the formulation of our toy model, there is 
no meaningful way for the mass-generating phase tran¬ 
sition to occur very much later than the critical over¬ 
damped/underdamped transition for either of our two 
fields: the lighter field is massless prior to the mass¬ 
generating phase transition and as such remains over¬ 
damped during this entire period, while the heavier field, 
either formally overdamped or underdamped, has no am¬ 
plitude of its own until the mass-generating phase tran¬ 
sition. For this reason, in this section we shall focus 
on those regions of parameter space in which our fields 
become underdamped near the time at which the phase 
transition begins. 

The most complete way of surveying physics within 
this regime is to vary both tq and Aq, with tq limited 
to the range tq ^ where each is implicitly de¬ 
fined (as in Fig. 2) as the time at which 3iL = 2Xi. This 
method, which by construction surveys all possibilities, 
thus involves variations in two independent parameters. 
However, another way of surveying many aspects of the 
physics within this regime is to fix tq, and instead to vary 
the width of the phase transition Aq- By considering 
fiducial values of tq which are carefully chosen with re¬ 
spect to the fnfj values, we can reach situations in which 

are near tq or just prior to it. These values of 
are illustrated in Fig. 9, and we see that this method also 
allows us to reach the desired values of It is there¬ 
fore this latter approach that we shall follow in the next 
two sections. 

Given our choices for tq and Aq, we can then cal¬ 
culate the late-time total energy density p for different 
values of the rhfj parameters — i.e., for different values 
of Wsum, Am^, and (f, 0max)- Our goal is to understand 
the effects on the total late-time energy densities p that 
emerge when a non-zero mixing, parametrized by and 
a non-zero width Aq for our mass-generating phase tran¬ 
sition are present simultaneously. 

Our results are shown in Figs. 10 through 12: In 
Figs. 10 and 11 we assume tg = 1 (for radiation- 



(j>\g in a matter-dominated universe, expressed as a number 
of widths Ag prior to the central phase-transition time tq 
and plotted as a function of Ag/tq for different values of tq. 
For tg ~ 1, the critical damping transition essentially occurs 
during the central portion of the mass-generating phase tran¬ 
sition for all Ag; for large tg, by contrast, the critical damp¬ 
ing transition occurs just prior to the mass-generating phase 
transition. These two cases thus survey our main regions of 
interest in this paper. 


and matter-dominated universes respectively), while in 
Fig. 12 we assume tq 1 (for which matter- and 
radiation-dominated universes yield the same results). 
In each of these figures we survey values of the late-time 
mixing saturation ^ within the range 0 < ^ < 1; note that 
we refrain from considering the actual limiting case with 
^ = 1 when discussing late-time quantities such as the 
late-time energy density p because the lighter field (fxg 
remains massless at late times when ^ = 1 (as evident 
from Fig. 6) and therefore never technically enters the 
asymptotic underdamped region which characterizes our 
definition of “late” times. We also note that although 
the left, center, and right panels within each figure cor¬ 
respond to different values of the maximum allowed mix¬ 
ing angle 0max, we see that it is only in terms of the 
mixing saturation ^ that we can make sensible compar¬ 
isons across the different panels. Indeed, although a fixed 
change in ^ within a panel with large 0niax corresponds 
to a much larger change in the absolute mixing angle 9 
than it does within a panel with small 0max, we see that 
it is only the changes in f rather than 6 which can be 
compared meaningfully across panels. For this reason, 
even though the value of 0max is important in order to 
express our results in terms of absolute mixing angles, 
in the following we shall simply describe our mixings as 
small or large depending on the corresponding values of 













13 



FIG. 10: Behavior of the total late-time energy density p(Ag,C) when mixing effects, parametrized by are combined with a 
non-zero width Aq for the mass-generating phase transition. For these plots we adopt a relatively small fiducial time tg = 1.0 
and assume a radiation-dominated universe; we also hold msum = 4 and take Am^ = —3, 0,-1-3 for the left, center, and right 
panels, respectively. In each case we plot the late-time energy density p as a function of the phase-transition width Aa for 
different values of the mixing saturation where in each panel p is normalized to its value for Aq = 0 and 5 = 0. 






FIG. 12: Same as Figs. 10 and 11, except that we now adopt a relatively large fiducial time tq = 10^. For such large values of 
TG, the results are nearly the same for both radiation- and matter-dominated universes. Note that unlike Figs. 10 and 11, the 
vertical axes here are logarithmic. 
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These figures illustrate the effects of turning on a fi¬ 
nite width Ag for the mass-generating phase transition in 
conjunction with non-zero mixing between our two fields. 
For the small -tq regime plotted in Figs. 10 and 11, we ob¬ 
serve that small mixing actually enhances the late-time 
energy density p. By contrast, we see that large mix¬ 
ing actually suppresses the late-time energy density p. 
Indeed, we see from Figs. 10 and 11 that these effects 
(in both directions) are more pronounced for matter- 
dominated universes than radiation-dominated universes 
and for situations in which Am? is positive rather than 
zero or negative. 

For the large-rc regime plotted in Fig. 12, by contrast, 
we see that matter- and radiation-dominated universes 
give rise to identical results. All late-time energy den¬ 
sities p, regardless of the mixing saturation share a 
common value when Ac = 0. However, unlike the small- 
T (3 case, the late-time energy densities experience no en¬ 
hancement at all, even in the absence of mixing = 0). 
Indeed, p experiences only a suppression for non-zero 
Ag — a suppression which, given the logarithmic p axes 
within Fig. 12, grows much more severe than it was for 
small tg- Moreover, we see that this suppression of the 
late-time energy density p is stronger for negative Am? 
(j.e., for negative fj) than for zero or positive — a feature 
which is completely reversed relative to the small-rg case! 
Finally, we observe the emergence of a non-monotonic 
“oscillatory” behavior for p as a function of Ag when tg 
is large and when the mixing saturation ^ grows close to 
1. This feature is most pronounced when Am^ is pos¬ 
itive, but exists for all An?. These oscillations will be 
discussed in Sect. V, and indicate a strong sensitivity of 
the suppression factor relative to even small variations in 
the phase-transition width Ag- 

All of these results illustrate the dramatic conse¬ 
quences that ensue when we consider the timescales 
associated with both our mass-generating phase tran¬ 
sition and the mixing it generates. For example, we 
see from the Am^ > 0 plot in Fig. 11 that mixing has 
no effect on the total late-time energy density p when 
the mass-generating phase transition is rapid {i.e., when 
Ag = 0). Thus, it is only the existence of a non-zero 
phase-transition timescale which allows the mixing to 
leave a non-trivial imprint at late times! Moreover, 
we see from these hgures that the enhancements and 
suppressions experienced by p are typically quite large, 
stretching from 20% or 30% in the case of enhancements 
all the way to many orders of magnitude in the case of 
suppressions! These effects can thus have dramatic im¬ 
plications for the relative size of the corresponding slice 
of the “cosmic pie” — i.e., for the overall composition of 
the total late-time energy budget of the universe. 

It is also instructive to understand those limits of the 
above results in which one or the other of our two vari¬ 
ables f and Ag is taken to zero. We begin by focus¬ 
ing on the effects that emerge solely due to the presence 
of a non-zero width Ag for the mass-generating phase 
transition {i.e., the effects that occur when ^ = 0). Re- 
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FIG. 13: Contours of the total late-time energy-density en¬ 
hancement/suppression factor p(Ag)/p(0), plotted in the 
(tg,o-) plane, assuming a matter-dominated universe with¬ 
out mixing = 0) and moo = 1- The thick black line running 
vertically and along the cr = 0 axis is the contour for which we 
have neither enhancement nor suppression. The correspond¬ 
ing contours for a radiation-dominated universe are similar. 
We observe that in the absence of mixing, a finite width for 
the mass-generating phase transition tends to modify the late¬ 
time energy density p(cr) compared to what it would have 
been for an instantaneous phase transition; this suppresses 
p((t) rather significantly for tg ^ 5, but enhances p(cr) for 
TG 1. Remarkably, these results are generally insensitive 
to (7 as long as cr > 0.1, and follow an approximate scaling 
behavior p ~ in this region. 


suits are shown in Fig. 13, where we plot contours of 
the total normalized late-time energy density p(o')/p(0) 
within the {tg,<t) plane, where u = Ag/(\A'T'g) as in 
Eq. (2.13). Generally, we see from Fig. 13 both suppres¬ 
sion and enhancements in the late-time energy density 
are possible, depending on whether tg ^ 5 or tg ^ 1 re¬ 
spectively. These results therefore explain the behaviors 
of the ^ = 0 curves in Figs. 10, 11, and 12. Furthermore, 
we see from Fig. 13 that p{Ag) is largely insensitive to 
a, and instead depends almost exclusively on the value of 
Tg- Indeed, for tg ^ 1, we see from Fig. 13 that p{Ag) 
scales approximately as Tq^. This in turn implies that 
p{Ag) Ag^ in this region. This observation will be 
discussed further in Sect. IV. 

Conversely, we can also study the effects that arise due 
to mixing alone {i.e., with Ag = 0). It turns out that 
when Ag = 0, the late-time energy density p can actu¬ 
ally be calculated analytically for arbitrary mixing; the 
general result is given in Eq. (A8). In Fig. 14, we choose 
values Aw? = 3 and mgum = 4, and plot the correspond¬ 
ing late-time energy density p{f) as a function of the 
late-time mixing saturation normalized to its unmixed 
value p(^ = 0), for a variety of different tg- For ^ = 0, 
this corresponds to a situation in which our two original 
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FIG. 14: Suppression of the late-time energy density p due to 
mixing effects alone, with Ag = 0, in a radiation-dominated 
universe. We assume an instantaneous phase transition at 
T = tq, and plot the corresponding late-time energy density 
p(5) as a function of the late-time mixing saturation nor¬ 
malized to its value at ^ = 0, for different choices of tq- We 
see that in all cases the late-time energy density p(^) ex¬ 
periences a suppression which grows increasingly severe as 
the mixing saturation is increased. As discussed in the text, 
the magnitude of this effect increases as p = Am ^—>■ 1 
(i.e., as fhoo —> 0), and is entirely absent if the phase transi¬ 
tion instead occurs during a matter-dominated epoch or more 
generally if re ^ 1. 


unmixed fields 0o,i remain uncoupled, with all of the en¬ 
ergy density remaining in ipo; by contrast, as ^ increases, 
increasing amounts of energy density are shared between 
our two fields during the mass-generating phase transi¬ 
tion. Surprisingly, we see from Fig. 14 that increasing 
^ results in an enhanced dissipation of the total energy 
density, so that the total late-time energy density p(^) 
is suppressed relative to what it would have been in the 
absence of mixing. Indeed, we see from Fig. 14 that this 
suppression is strongest for relatively small tq and rela¬ 
tively large mixing saturations Moreover, as evident 
from the calculation in Appendix A, this effect exists in 
all cases except for a matter-dominated epoch. This ef¬ 
fect was first observed in Ref. [35] within the context of an 
infinite tower of Kaluza-Klein axion modes, where it was 
exploited in order to permit a higher-dimensional loosen¬ 
ing of the usual four-dimensional overclosure bounds on 
the Peccei-Quinn scale /pq. However, we now see that 
this effect is completely general, and persists even when 
only two modes are involved. 

We thus conclude that for phase transitions occurring 
during a radiation-dominated epoch, mixing alone in¬ 
duces a significant suppression in the late-time energy 
density. Indeed, this is true even if the phase transition is 



FIG. 15: Suppression of the late-time energy density p due 
to mixing effects alone, with Ag = 0, for arbitrary universes 
parametrized by k. As in Fig. 14, we assume an instantaneous 
phase transition at r = tg, but we now plot the corresponding 
late-time energy density p{f) as a function of n, normalized 
to its value at ^ = 0, for different choices of We see that 
mixing causes a suppression of the late-time energy density 
within K < 2 universes but an enhancement within k > 2 uni¬ 
verses — indeed, these effects disappear only in the special 
case of a matter-dominated universe. These effects can there¬ 
fore produce significant modifications of the final late-time 
energy density, even for instantaneous mass-generating phase 
transitions. 


treated as instantaneous. This effect grows in significance 
as Amf —>■ and ^ t 1. This effect is also largest 

(i) (0 

when Tq ^ , and ultimately vanishes for tq . 

Note that all of these observations are consistent 
with the plots shown in Figs. 10, 11, and 12. For 
Ac = 0, we see that all of the plotted curves begin at 
p(A(3,^)/p(0, 0) = 1 except for those in Fig. 10, where 
the mere fact of having ^ 0 induces a suppression of 

the late-time energy density, even for Aq = 0. 

Finally, along the same lines, it is interesting to con¬ 
template what happens for more general universes be¬ 
yond those that are radiation-dominated. The corre¬ 
sponding results for universes with general values of n 
are shown in Fig. 15. We see that mixing alone indeed 
produces a suppression of the late-time energy density for 
universes with k < 2, but this suppression actually be¬ 
comes an enhancement for universes with k > 2! More¬ 
over, as expected, both effects become stronger as the 
mixing saturation f increases — strong enough to change 
the late-time energy density by factors of two or three or 
even more. Indeed, it is only for a matter-dominated 
universe (corresponding to k = 2) that this effect disap¬ 
pears. 
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IV. INDIVIDUAL LATE-TIME ENERGY 
DENSITIES 

In the previous section, we focused on only one quan¬ 
tity: the total late-time energy density p- However, an¬ 
other important feature is the actual distribution of this 
total energy density p amongst the two fields of our sys¬ 
tem. We shall now proceed to study this issue. 

In Fig. 16, we show the behavior of the individual late¬ 
time energy densities associated with (the lighter 
field) and (j)x^ (the heavier field) for a mass-generating 
phase transition occurring at tq = 10"^. It is therefore 
the sum of these two energy densities which produces 
the results in Fig. 12; recall that these results are the 
same for matter- and radiation-dominated universes. It 
turns out that there are many features illustrated within 
Fig. 16 which will be important for our future results. 
We shall therefore step through these features, one by 
one. 

We begin by concentrating on the special case without 
mixing (he., with ^ = 0). This will allow us to study the 
effects of the non-zero phase-transition width Aq for a 
single component alone without the complications due to 
mixing. For Am^ > 0, all of the energy density arising 
due to the mass-generating phase transition accrues to 
the lighter field (j)Xg- For Am^ < 0, by contrast, all of 
the late-time energy density is associated with the heavier 
field <j)Xi, as discussed at the end of Sect. II. 

In either case, we observe from Fig. 16 that the cor¬ 
responding late-time energy density px remains essen¬ 
tially constant (he., independent of A^j) for A^ < 2'k/X. 
This makes sense, as the effects of introducing a non-zero 
width for our phase transition will be essentially invisi¬ 
ble if the mass generation occurs more rapidly than the 
natural timescale of oscillations of our (underdamped) 
field 4>x- In other words, for Aq ^ Stt/A, the process 
of pumping energy density into our system occurs with 
what may be considered to be maximum efficiency, since 
the phase transition appears to be effectively instanta¬ 
neous with respect to the natural oscillation timescale. 
However, for Aq ^ Stt/A, the field oscillations tend to 
compete against the process of mass generation. As a re¬ 
sult, the mass-generating phase transition is less efficient 
in pumping energy density into our system, thereby in¬ 
ducing a suppression in the late-time energy px which 
ultimately scales as an inverse power of Ac- 

Px{Ag,0) ~ 1/Ag for Ag » 27 r/A . (4.1) 

Thus, once Ag exceeds the natural oscillation period of 
our mass eigenstate, increasing the width of the phase 
transition has the effect of introducing a power-law sup¬ 
pression of the corresponding late-time energy density. 

These results hold for the unmixed scenarios with 
^ = 0. However, for the lighter fields (j)Xo (correspond¬ 
ing to the top row of Fig. 16), we see that the above 
asymptotic behavior continues to hold regardless of the 
value of the mixing: 

Pao(^G>?) ~ 1/^G for Ag » 27r/Ao . (4.2) 


Thus, we see that increasing the width of the phase tran¬ 
sition continues to be associated with a power-law sup¬ 
pression of the late-time energy density of the lighter field 
— even in the presence of non-zero mixing. 

Turning on a mixing between our two fields also has 
a number of other important effects on their individual 
late-time energy densities. For the lighter fields (jixg (as 
considered along the top row of Fig. 16), the most promi¬ 
nent effect is of course the set of very strong oscillations 
which are induced for very large mixing saturations ^ ^ 1 
and “intermediate” widths 1 ^ Ag ^ 100. We shall defer 
our discussion of these oscillations until the next section, 
but we see from Fig. 16 that these oscillations are rel¬ 
atively large, occasionally enhancing the corresponding 
late-time energy densities above what they would have 
been in the absence of mixing by several orders of mag- 
nitudel Indeed, this enhancement of the corresponding 
late-time energy densities above their ^ = 0 values per¬ 
sists even for values of Ag which lie beyond the actual 
oscillations themselves. 

As a result, we see that in the presence of a non-zero 
width for the phase transition, mixing has a general ten¬ 
dency to enhance the late-time energy density of the 
lighter field — all without disturbing the power-law sup¬ 
pression discussed above. Indeed, for the lighter fields 
(fxg , the effects of the non-zero transition width Ag and 
the non-zero mixing saturation ^ tend to pull in opposite 
directions: the width tends to suppress the correspond¬ 
ing late-time energy density while the mixing tends to 
enhance it. As evident from the plots along the top row 
of Fig. 16, the nature of the net result (either an overall 
enhancement or overall suppression) therefore depends 
non-trivially on the precise values of Ag and f involved. 

By contrast, for the heavier fields 4>Xi (for which the 
corresponding energy densities are plotted along the bot¬ 
tom row of Fig. 16), the effects of non-zero mixing are 
quite the opposite: the enhancement discussed above is 
gone, and instead there is now an additional suppression 
which helps to drive the corresponding energy densities 
to even smaller values, as functions of Ag, than we had 
for the lighter fields! 

It is important to understand the origins of this ad¬ 
ditional suppression. Unlike the lighter field (fxo^ which 
obtains its energy density directly from the phase tran¬ 
sition through the generation of a non-zero mass, the 
heavier field (fxi actually gains its energy density when 
it gains an overall field amplitude through mixing with 
the lighter field. The magnitude of the resulting energy 
density is then governed by two factors: the magnitude 
of the field amplitude generated, and the masses that are 
also generated by the phase transition. 

Thus, while the natural timescale governing the rate at 
which the energy density is originally pumped into the 
lighter field is given directly by Ag — and thus by the 
function h(r) plotted in Fig. 3 — the natural timescale 
governing the rate at which the energy density is pumped 
into the heavier field is governed not only by Ag but also 
by the rate of change of the mixing angle 9 plotted in 
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FIG. 16: Behavior of the individual late-time energy densities ^ {Ao,0 when mixing effects, parametrized by are combined 
with a non-zero width Ag for the mass-generating phase transition. For these plots we adopt the large fiducial time tq = 10"^ 
utilized for Fig. 12, which leads to identical results for radiation- and matter-dominated universes. The top row shows results 
for P)^^ (associated with the lighter state in all cases), while the bottom row shows results for p^^ (associated with the heavier 
state); we have also held mLm = 4 and take Am^ = —3,0,-1-3 for the left, center, and right columns, respectively. In each 
case we plot the late-time energy density p^ as a function of the phase-transition width Aq for different values of the mixing 
saturation normalized to its value for Ac = 0. Note that we refrain from plotting curves for ^ = 0 (and instead plot curves 
for ^ = 0.05) in cases where the corresponding energy densities vanish. 


Fig. 8. In analogy with Eq. (2.12), we may even define 
a corresponding width Ag for the mixing angle 9 via the 
slope of 9{t) at its midpoint 9/2: 


Ag 


(9/9) 


T6 


(4.3) 


here rg, as defined in Eq. (2.33), is defined analogously 
to tg as the time at which 9 = 9/2. A rough analytical 
approximation for Ag when fj < 0 is given by 


Ag 



for T] < 0 . 


(4.4) 


The magnitudes of Ac and Ag are compared in Fig. 17, 
from which we see that Ag never exceeds Aq- Indeed, 
we see from Fig. 17 that Ag can be much smaller than 
Ag for 7) < 0 and ^ ^ 1. We thus see that during our 
mass-generating phase transition, the mixing angle 9 al¬ 
ways changes at least as rapidly as do the eigenvalues A^, 
and can in fact under certain circumstances change even 


more rapidly, the latter despite the fact that both sets of 
changes arise due to the same phase transition. 

This observation is directly relevant for the manner in 
which the heavier field accrues its energy density from the 
mass-generating phase transition. As such, this feature 
is thereby directly relevant for the resulting behavior of 
the late-time energy density . There are three distinct 
cases to consider. 

• Ag 2'k/\\\ In this case, the oscillations of the 
heavier field are so slow during the phase transi¬ 
tion that all aspects of the phase transition appear 
effectively instantaneous. This is true for the gener¬ 
ation of the amplitude of the heavier field as well as 
the change in its mass. Energy density is thereby 
delivered to the heavier field with maximum effi¬ 
ciency, in a manner which is completely insensitive 
to the non-zero timescales Ag and Ag associated 
with the phase transition and which is thereby pro¬ 
tected against all of their associated dissipating ef¬ 
fects. This behavior is evident for sufficiently small 
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FIG. 17: The ratio of the two timescales Ag and Ag, plotted 
as a function of fj for different mixing saturations In all 
cases, we see that Ag < Ag, with Ag <C Ag for fj < 0 and 
^ <C 1. The dots represent the approximate analytical result 
in Eq. (4.4). The distinction between Ag and Ag has impor¬ 
tant implications for the energy density associated with the 
heavier field <)>Ai , since Ag governs the rate at which the mass 
of this field is modified during the phase transition while Ag 
governs the rate at which the amplitude of this field is gener¬ 
ated as a result of its mixing with the lighter field. 


Ag within the plots along the lower row of Fig. 16. 

• Ag < 27r/Ai < Ag- In this case, the oscillations of 
the heavier field are sufficiently slow that the gener¬ 
ation of its field amplitude appears to be effectively 
instantaneous; these oscillations are nevertheless 
sufficiently rapid that full oscillations can occur 
during the change in field mass. The net result is 
that the heavy-held energy density experiences 
the same suppression as does the light-held energy 
density Paq’ with both depending non-trivially on 
Ag and ultimately scaling inversely with Ag. Note 
that this case exists only for situations in which Ag 
is signihcantly smaller than Ag, which according 
to Fig. 17 tend to emerge only for Am^ < 0 and 
relatively small mixing. 

• 27r/Ai < Ag: In this case, the held oscillations of 
the heavier held are sufficiently rapid that full os¬ 
cillations are occurring not only during the change 
in its mass but also during the generation of its am¬ 
plitude. This latter feature leads to an additional 
source of suppression for the late-time energy den¬ 
sity beyond what the lighter held experiences, and 
causes to exhibit an even more dramatic sup¬ 
pression as a function of Ag than that exhibited by 
PAq ■ This behavior is clearly evident for sufficiently 
large Ag in the plots along the lower row of Fig. 16. 


The above discussion describes the behaviors of the in¬ 
dividual components Pa as functions of Ag and f. How¬ 
ever, our hnal task is to determine how much each of 
these individual components contributes to the total en¬ 
ergy density p. Of course, to do so requires that we under¬ 
stand not just the intrinsic behaviors of these individual 
components as functions of Ag and but also the rela¬ 
tive sizes of these individual components. In other words, 
we need to understand the relative normalization of the 
curves for Paq in ths top row of Fig. 16 relative to those 
for Pai in the bottom row. However, this relative normal¬ 
ization may easily be determined in the Ag = 0 limit, for 
which analytical results are given in Appendix A. Indeed, 
use of the virial approximation A^rg ^ 1, which is valid 
for the plots in Fig. 16, yields the ratio 

Pao(O)C) ^ Apcos^g ^ 2 A -rjsec2e \ 

mo " A^cos2 0 + A^sin^0""°' V 1-^ ) 

(4.5) 

as the fraction of the total late-time energy density re¬ 
maining in the lighter field. As expected, this quan¬ 
tity evolves monotonically from unity at 0 = 0 to zero 
at 0 = 0max [where 0max is defined in Eq. (2.26)]. It 
may thus be greater or less than 50%, depending on the 
eigenvalues and mixing angles in question. Note that in 
the limit of small mixing (which corresponds to 0 « 0 
or 0 « 7r/2), the result in Eq. (4.5) implies that whether 
the bulk of the energy density winds up at late times 
associated with the lighter or heavier field depends sim¬ 
ply on the sign of p rather than its magnitude. (Re¬ 
call, in this connection, that p —>■ —p implicitly changes 
0 ^ 7r/2- 0.) 

The result in Eq. (4.5) holds only for Ag = 0. How¬ 
ever, as Ag grows larger, this ratio will change. For the 
total energy densities plotted in Fig. 12, the correspond¬ 
ing fractional results are shown in Fig. 18 as a function 
of Ag. It is easy to understand the behavior shown in 
Fig. 18. Regardless of the fraction of the total late-time 
energy density associated with the lighter field ^Aq when 
Ag = 0, we have already seen that the energy densities 
associated with the heavier fields are typically more sup¬ 
pressed as a function of Ag than are the energy densi¬ 
ties associated with the lighter fields. Thus, regardless of 
the late-time energy-density configuration when the phase 
transition is instantaneous, we see that increasing the 
width of the phase transition has the effect of throwing 
an increasingly large share of the total late-time energy 
density into the lighter field. Indeed, in some cases, we 
see that we can entirely reverse the distribution of the 
total energy density from the heavier field to the lighter 
field, simply by adjusting the timescale over which the 
phase transition occurs! 
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FIG. 18: The fraction of the total late-time energy density p{Ag,^) which is associated with the lighter field 4>\q, plotted as 
functions of Ag for different values of with tq = 10"^. Regardless of the small-A g values of this fraction [as given analytically 
in Eq. (4.5)], we see that increasing the width Ag of the phase transition generally has the net effect of transferring more and 
more of the total energy density into the lighter field. Thus, simply by adjusting the width of the phase transition, we see that 
we can often entirely reverse the distribution of the late-time energy density. 


V. PARAMETRIC RESONANCE 

We now turn to what is perhaps the most prominent 
feature within the plots shown in Figs. 12, 16, and 18: 
the appearance of non-monotonicities in the late-time en¬ 
ergy densities p\ associated with the lighter fields (j)\g 
as functions of the phase-transition width Aq. As in¬ 
dicated in Sect. Ill, these oscillations grow particularly 
large when tq is large and when the mixing saturation 
^ grows close to 1. These non-monotonicities indicate 
that within this region of parameter space there exists an 
extremely strong sensitivity of the late-time energy den¬ 
sity Pao fo even small variations in the phase-transition 
width Ag. Moreover, as evident from Fig. 16, these non¬ 
monotonicities can often enhance Pao by several orders 
of magnitude! 

These non-monotonicities are yet another consequence 
of the interplay between the width of our phase tran¬ 
sition and the mixing it induces. It is not difficult to 
understand the origin of these oscillations: ultimately 
they are parametric resonances triggered by our mass¬ 
generating phase transition. Recall that a parametric 
resonance generically occurs when the mass of an oscilla¬ 
tor itself exhibits an oscillatory behavior whose frequency 
is approximately twice the natural frequency of the oscil¬ 
lator. (An example of this is a child “pumping” a swing 
by alternatively standing and squatting on the swing seat 
at the correct moments during the period of the swing.) 
In such systems, the amplitude of oscillation grows expo¬ 
nentially and without bound unless the oscillator also has 
a frictional damping term. A similar behavior can also 
emerge if the friction term (rather than the mass term) 
experiences the oscillatory behavior. 

At first glance, it may seem that our system of coupled 
scalar fields has no means of experiencing a parametric 
resonance. The equations of motion of our model are 
given in Eq. (2.7), and although both the Hubble damp¬ 


ing term H{t) and the mass matrix in our model 
are time-dependent, the Hubble term is monotonically 
falling. Likewise, the /i(r)-function — which governs the 
time-development of our mass matrix, as in Eq. (2.9) — 
is also monotonic. 

However, this is where the mixing plays a critical role. 
Because of the non-zero mixing inherent in this sytem, 
the aspects of the mass matrix which are important are 
not its individual components but rather its eigen¬ 
values Xi- Indeed, even though the mass-matrix compo¬ 
nents mfj{T) are all monotonic as functions of time r, the 
eigenvalues Xi need not be. These eigenvalues were plot¬ 
ted as functions of time in Fig. 7, and we see that for the 
smaller eigenvalue Aq the mixing induces a level-repulsion 
that often results in a non-monotonic “pulse”. This pulse 
is more prominent for < 0 than for Am^ > 0, but it 
is ultimately a generic feature of the behavior of the lower 
eigenvalue for sufficiently saturated mixing. Indeed, this 
pulse grows increasingly dramatic as ^ I, and has a 
width governed by Aq. In general, it is straightforward 
to show that a pulse will always appear for the lightest 
eigenvalue if the extra contribution to the squared-mass 
matrix that comes from the mass-generating phase tran¬ 
sition has negative determinant — i.e., detmT < Q. 

Strictly speaking, a single pulse does not oscillatory 
behavior make. However, during the relevant time inter¬ 
val near tq, this pulse may be regarded as one oscillation 
within a full sinusoidal pattern. This situation is illus¬ 
trated schematically in Fig. 19. Of course, variations in 
Ac directly affect the width of this pulse, and thereby 
change the effective pulse frequency. As a result, there 
will exist a single value of Aq for which this frequency 
is exactly twice the natural frequency of our oscillator, 
and for which our primary parametric resonance emerges. 
Even greater values of Aq then correspond to the higher 
harmonics of this resonance. 

This, then, is the origin of the parametric oscillations 
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FIG. 19: A schematic illustration of the “pulse” (solid blue 
curve) experienced by the smaller eigenvalue Ao(t) near 
T ~ Tp and a full sinusoidal function (dashed red curve) to 
which it may be approximated in the vicinity of the pulse. 
As discussed in the text, this approximation allows us to ex¬ 
tract the effective frequency uJes associated with the pulse and 
thereby determine the mathematical condition under which 
an n*’'-order parametric resonance occurs. 


apparent in the plots shown in Figs. 12, 16, and 18. 
As such, this entire phenomenon is a prime example of 
the interplay between the mixing ^ and the width of 
the mass-generating phase transition: the non-zero mix¬ 
ing produces the level repulsion that leads to the pulse, 
while the non-zero width of the phase transition endows 
this pulse with the specific width/frequency needed for 
it to potentially trigger an actual resonance. Moreover, 
with this explanation, we now understand why these res¬ 
onances occur only for certain values of Ag, and only 
for the lighter field. This also explains why these reso¬ 
nances grow stronger as ^ ^ 1. Furthermore, although 
this pulse exists for both < 0 and Am^ > 0 (as ev¬ 
ident in Fig. 7), we have already learned in the previous 
section that for the intermediate values of Ag relevant 
for the parametric resonance, the energy density of the 
lighter field contributes a greater fraction of the total 
energy density for Ato^ > 0 than it does for Am^ < 0 . 
This behavior is illustrated, for example, in Fig. 18. It 
is for this reason that the total energy density p exhibits 
a stronger parametric resonance for Am^ > 0 than for 
Am^ < 0 (as evident in Fig. 12), even though the in¬ 
dividual component ]}xg exhibits a stronger parametric 
resonance for Am^ < 0 than for Am^ > 0 (as evident in 
Fig. 16). 

In general, it is easy to determine the values of Ag for 
which such parametric resonances occur. As indicated 
above, we imagine that Ag(T) experiences a “pulse” cen¬ 
tered at r = Tp, as sketched in Fig. 19. Within the region 
of the pulse, we then approximate the behavior of the 
pulse as part of a sinusoidal function with an effective 


frequency Wes: 


for r 



~ 5Ao(Tp) |l -b COs[Weff(T - Tp)] 

From Eq. (5.1) we then find that 



(5.1) 


(5.2) 


A pulse with this frequency will lead to an n*^-order 
parametric resonance only if LOeS is ( 2 /n) times the nat¬ 
ural oscillation frequency of our system. The case with 
n = 1 produces the primary parametric resonance, while 
the cases with n S Z > 1 produce its higher-order har¬ 
monics. However, near r « Tp, the natural oscillation 
frequency associated with the lightest held is nothing but 
Ao(Tp), since it is the mass of the held which drives its 
oscillations. We thus must demand 


^eff 


2 

n 


■ >^0{Tp) ■ 


(5.3) 


In other words, a parametric resonance will occur only 
when two timescales are balanced: the timescale l/weff 
of the pulse, and the timescale l/Ao(rp) of the held oscil¬ 
lations in the vicinity of the pulse. Combining this with 
Eq. (5.2), we then obtain a condition for an n*^-order 
parametric resonance: 





n e Z+ . (5.4) 


Variations in Aq will modify the value of the left side of 
this equation. Thus, there exist a discrete set of values 
Aq'^ for which this equation can be satished, one for each 
value of n. These are then the phase-transition widths 
for which parametric resonances exist. 

The solutions for Aq'^ can be obtained numerically, 
and our results are shown in Fig. 20. We see that the 
values of A[?^ generally increase with n, as expected. 
Likewise, these critical widths also increase as functions 
of fj and ultimately diverge as 77 —?> 1. Moreover, super¬ 
imposed on this plot are contours showing the values of 
27r/Ai in relation to Ag. As noted in Sect. IV, the smaller 
27r/Ai is in relation to Ag, the more suppressed is the 
energy-density contribution of the heavier field relative 
to that of the lighter held and hence the more the lighter 
held (which feels the parametric resonance) dominates 
the total energy density. This then conhrms our previ¬ 
ous expectation that the parametric resonance, through 
stronger within when Am^ < 0 , is nevertheless more 
pronounced within p when Ato^ > 0 . 

Thus, we conclude that parametric resonances not only 
occur, but can also dominate the total energy density of 
our system. Such parametric resonances can distort the 
resulting late-time energy densities by several orders of 
magnitude, and thus may play an important role in early- 
universe cosmology. 
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FIG. 20: The phase-transition widths (n=l,2,3) at 
which the first three parametric resonances occur, plotted as 
functions oirj. Contours showing 27r/Ai in relation to Ae are 
also superimposed. The region with large fj is thus one in 
which our parametric resonances occur while Ag ^ 27r/Ai — 
i.e., a region in which the total late-time energy density p is 
dominated by the contribution from the lighter field experi¬ 
encing the parametric resonance. 


VI. RE-OVERDAMPING 

We now turn to another novel feature which results 
from the confluence between mixing and a finite phase- 
transition width: a phenomenon which we refer to as 
“re-overdamping”. 

As discussed at the beginning of Sect. II, any scalar 
field (j) with mass X{t) will experience overdamped behav¬ 
ior (and thus function as vacuum energy) at times t for 
which 3i?(<)/2 > X{t). By contrast, such a field will expe¬ 
rience underdamped behavior (and thus function as mat¬ 
ter) when 3iJ(t)/2 < X{t). Of course, H{t) ^ 1/t in any 
post-inflationary epoch. Thus, in situations for which 
X(t) is either constant or monotonically increasing, as il¬ 
lustrated in Fig. 2, there exists a single, relatively short 
era of times t for which 3iF(t)/2 « A(t). Indeed, this 
is the period, as discussed in Sect. II, during which our 
scalar held (p is transitioning from overdamped to under¬ 
damped behavior. After this transition, our held then 
remains underdamped for all future times. 

This is the standard situation, and assumes that X{t) is 
either constant or monotonically increasing. However, as 
we discussed in Sect. V, this need not be true when multi¬ 
ple helds mix in the presence of a mass-generating phase 
transition of non-zero width. Indeed, as we have seen in 
Sect. V for our two-component toy model, the mass of 
our lighter held can experience a “pulse”, as illustrated 
in Fig. 19; this is the phenomenon underlying the para¬ 


metric resonances discussed in Sect. V. In the presence of 
a pulse, however, the 2>H{t)/2 curve may cross the X{t) 
curve not just once but three times, resulting in a held 
which begins overdamped, then becomes underdamped, 
but then experiences a re-overdamped phase before even¬ 
tually settling back into a hnal underdamped state. This 
behavior is illustrated in Fig. 21. Moreover, as sketched 
in Fig. 21, the period during which the held exists in a 
re-overdamped phase can be of extremely long duration 
— even relative to the width of the pulse that produced 
it. Indeed, as we shall discuss below, this duration can 
become arbitrarily long, depending on the values of the 
underlying input parameters of our model. As a result, 
the re-overdamped phase can persist throughout a sizable 
portion of cosmological history. 

It is important to stress that unlike the original over¬ 
damped phase, this re-overdamped phase of the theory 
should not be associated with vacuum energy. Indeed, at 
the moment of transition into the re-overdamped phase, 
oscillatory behavior ceases but our field value re-enters 
the overdamped phase carrying a non-zero velocity 4>. 
Thus, the re-overdamped phase is one in which our field 
value has a non-zero but steadily decreasing first deriva¬ 
tive. In this sense, then, the re-overdamped phase rep¬ 
resents a truly different behavior for our held values — 
one which is neither vacuum energy (constant held val¬ 
ues) nor matter (oscillatory held values), but something 



FIG. 21: A schematic illustrating the behavior of the lower- 
mass eigenvalue Xo{t) in the presence of a “pulse”, plotted as a 
function of time with the Hubble parameter 3H(t )/2 superim¬ 
posed. The corresponding field 4>\g therefore experiences not 
just the usual transition from overdamped to underdamped 
behavior, as in Figs. 1 and 2, but also a novel re-overdamping 
transition. As shown, this re-overdamped phase can persist 
for considerably longer than the width Ag of the pulse that 
produced it, and thus can have a significant impact on the 
resulting cosmological history. 
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FIG. 22: The lighter field (j)\g (top row) and corresponding total energy density p (bottom row), plotted as functions of time 
for situations (left to right) in which ^ 1 (or e = 1 — ^ 0). Note that the fields and corresponding energy densities are 

normalized so that (j>\g = 1 for t tq. Time intervals during which the system is overdamped are shaded (blue), while the 
width Ag of the original mass-generating phase transition is also indicated (red). As the mixing saturation increases towards 
its maximum value (a point of enhanced symmetry for the system), we see that a re-overdamped phase emerges with increasing 
duration. For comparison purposes, superimposed on each plot are the field values and/or total energy densities that our system 
would have exhibited in the Ag —>■ 0 limit (green); recall that in this limit it is the energy density of the heavier field which 
dominates the total energy density. We thus see that the re-overdamped phase can easily persist throughout a large period of 
cosmological history, with a phenomenology that is quite different from those of the traditional overdamped or underdamped 
phases. 


completely different which transcends our traditional ex¬ 
pectations for scalar fields in an expanding universe. 

The re-overdamped phase tends to emerge in systems 
in which the mixing is close to its maximum allowed 
value, with ^ 1 (or e = 1 — ^ 0). The resulting be¬ 

havior for field values and energy densities in this limit 
are illustrated in Fig. 22. As expected, we see that the re- 
overdamped phase has a phenomenology which is quite 
distinct from those of the traditional overdamped or un¬ 
derdamped phases. 

On the one hand, this re-overdamped phase can be con¬ 
sidered to be a mere “transient”. After all, it cannot last 
forever: as long as the field has non-zero mass at late 
times, this phase must always give way to a final under¬ 
damped phase. In other words, the re-overdamped phase 
can never represent the eventual end-state for our system 
as < —>■ oo. Moreover, even within this phase, the non¬ 
zero initial field velocity which characterizes its unique 
properties always eventually dissipates to ^ « 0 as long 


as this phase has sufficient duration. In such cases, this 
phase eventually comes to resemble vacuum energy. On 
the other hand, however, this attitude can be viewed as 
unnecessarily stringent. Even though the re-overdamped 
phase must always eventually give way to one which is 
underdamped, we have seen that it can, depending on 
the specific model under study, be made to persist for 
arbitrarily long duration. Likewise, this phase can begin 
with arbitrarily large initial velocity (j). Thus, in purely 
practical terms, this phase can be considered alongside 
“vacuum energy” and “matter” as an equally valid new 
behavioral phase that a scalar field can experience during 
a significant portion of cosmological evolution. 

As a final comment, we remark that this re¬ 
overdamping phenomenon also leads to a resonance-like 
behavior wherein quantities such as the late-time energy 
density experience a non-monotonic sensitivity to quanti¬ 
ties such as Aq- This is because the behavior of our field 
(j) during the re-overdamped epoch depends very sensi- 
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tively on the particular field velocity (j) that happened 
to exist at the precise moment the re-over damping tran¬ 
sition occurs. This gives rise to a periodic sensitivity 
of the final energy density to any quantity which affects 
the positioning of the re-overdamping transition moment. 
However, it turns out that the mathematical conditions 
underlying this periodicity are essentially the same as 
those underlying the parametric resonance discussed in 
Sect. V; this makes sense as both concern the fitting 
of (half-)integral numbers of field oscillations within the 
same eigenvalue pulse. Moreover, it turns out that the 
parametric resonance necessarily exists for all situations 
in which re-overdamping exists (even though the converse 
is not true). As a result, in such situations both features 
are inherent in the solutions to our equations of motion 
and the effects of this re-overdamping-induced oscillation 
are difficult to disentangle from those of the parametric 
resonance. 


VII. PHENOMENOLOGICAL EXAMPLE: 
A SECOND AXION 


As a phenomenological example of some of the main 
results of this paper, let us consider the effects that might 
come from adding a second axion to the standard QCD 
axion theory. As we shall see, the introduction of a 
second possible axion can severely distort the physics 
normally associated with the ordinary QCD axion, even 
if the second axion is itself associated with a relatively 
heavy Peccei-Quinn scale or is only mildly coupled to the 
first. 

Our theoretical set up is as follows. We assume that (j)o 
and (f>i are both axion-like particles — i.e., the pseudo- 
Nambu-Goldstone bosons associated with the breaking of 
global symmetries t/(l)o and C/(l)i at energy scales /o 
and fi , respectively. We assume that both fields are cou¬ 
pled to the QCD gauge fields via the chiral anomaly, with 
0(1) model-dependent coefficients and ^i. Moreover, 
in keeping with the assumptions of our toy model, we 
assume that at early times (po is massless while (pi has 
a non-zero mass M which possibly arises due to some 
earlier dynamics. Eventually, non-perturbative QCD in- 
stanton effects give rise to a temperature-dependent ef¬ 
fective potential V{pi,p 2 ) which takes the form 


V{T) 


Ma ^qcd(^) 


1 — cos 


f ^oP’o 
\ fo 


h 



where AQQj3/(327r^), where /iqcd( 7’) describes 

the temperature-dependence of this effective potential as 
it turns on (as illustrated in Fig. 4), and where 0 is 
the QCD theta-angle. Thus we see that the instanton- 
induced QCD confining phase transition plays the role of 
our mass-generating phase transition: it occurs at time 
tG = tQCB with width Aq = Aqcd = and 

gives rise to an additional contribution to the squared- 
mass matrix of our axion-like fields which at late times 


takes the form 


ml 


d'^V 


d(pid(pj 


( 0 ) 


hfj 


(7.2) 


where we have set = 1 for simplicity. 

Given this setup, we thus have two mass eigenstates, 
p\g and with late-time masses Aq and Ai respec¬ 
tively. As we survey the range from M = 0 to M = oo, 
we find that these late-time masses increase monotoni- 
cally within the ranges 


0 < Ao < Ma//o (7.3) 


and 


T^Ayy^ + jA < Ai < oo . (7.4) 

For reasons to be discussed shortly, we shall identify the 
lighter mass eigenstate as our QCD axion, which means 
that we shall identify Aq as the axion mass ma- We shall 
also define x = Ao/oZ/r^; thus 0 < a: < 1. Note that in 
this section we are keeping M explicit in our expressions; 
thus, all quantities such as Ai have their true physical 
mass dimensions. Our goal is then to understand how 
the late-time energy density associated with this field is 
influenced by the presence of the additional axion. 

Note that there are two physically distinct ways in 
which to realize a limit in which the second axion decou¬ 
ples from the first: either we can take M —>■ oo and leave 
fi /fo arbitrary, or we can take fi/fo —>■ oo and leave M 
arbitrary. In the first case, we find that x = 1 regardless 
of the value of fi. In the second case, by contrast, we 
find that x = 1 for all M > but x = Mfo/n\ < 1 

for all M < fJ.\/fo- Of course both decoupling limits 
smoothly merge together in the region for which both 
M and /i//o are taken to infinity. 

In order to quickly survey the expected phenomenolo- 
gies of this two-axion system as functions of {M, /o, /i} 
— and also to understand the properties of these differ¬ 
ent decoupling limits — let us calculate the correspond¬ 
ing values of fj and In terms of {M, fo, fi}, these are 
given by 

_ 1 + {fi/fo? [2{Mfo/^il)^ - 1] 

1 + ifi/for [mfo/^^lr + 1 ] ^^ 

and 


v/1 + 2(/i//o)2(M/o/Mi)2 ■ 

These quantities are shown in Fig. 23. Although both de¬ 
coupling limits result in ^ 0 (the absence of mixing), as 

expected, we immediately see that the decoupling limit 
with fi ^ fo and arbitrary M generally gives rise to all 
possible values of fj in the range — 1 < ?? < 1. By contrast, 
the decoupling limit with M —> oo gives rise to 77 = 1 re¬ 
gardless of the value of fi/fo- 
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FIG. 23: The quantities rj (left panel) and ^ (right panel), plotted as functions of M for different values of /i//o. Superimposed 
on each plot are contours of x (dashed black lines). The decoupling regions in each case correspond to taking /i//o —> oo 
for any M as well as taking M ^ oo for any /i//o; these smoothly connect to form the curves labeled “decoupling limit” on 
each plot. Note that we have ^ ~ 0 along the entire decoupling-limit line, as expected since the decoupling limit is one in 
which there is no mixing between the two axions. By contrast, the decoupling-limit line surveys all possibilities of fj. We also 
observe from either panel that x = 1 along the decoupling line for all M > fo, while x < 1 along the decoupling line for all 
M < Pa//o; this is easier to see in the right panel but is true as a general statement about the relationship between M and x 
in the decoupling limit. Thus, an infinite M implies that x = 1 but x = 1 implies only that M > /taZ/o- 


This is an important distinction because we have al¬ 
ready seen in Sect. IV that fj > 0 is the regime in which 
the majority of the energy density is found in the lighter 
field at late times, whereas ?) < 0 is the regime in which 
the majority of the energy density is transferred to the 
heavier field at late times. This behavior was in fact dis¬ 
cussed explicitly below Eq. (4.5) in the limit of an instan¬ 
taneous phase transition. Thus, it is only the decoupling 
limit with M —> oo and arbitrary /i//o for which the en¬ 
ergy density actually remains in the lighter held where 
it started, as we would expect for a proper decoupling 
limit. 

We shall therefore adopt the limit with M —>■ oo and 
arbitrary /i//o as our benchmark against which to mea¬ 
sure the effects that come from the presence of the second 
axion. Note that in the M —>■ oo limit, the physics of our 
system actually becomes independent of /i//o; we have 
already seen this behavior for the specihc quantities fj and 
^ in Fig. 23, but this feature indeed holds more generally 
for all properties of the system. 

Although we have been considering the lighter mass 
eigenstate (f>\g as our QCD axion, it is important to re¬ 
alize that the entire notion of a “QCD axion” no longer 
strictly applies. For example, there is no mass eigen¬ 
state in our model whose mass is inversely proportional 
to /o or /i, as might be taken to characterize a QCD ax¬ 
ion. Likewise, there is no mass eigenstate in this model 
which solves the strong CP problem; in particular, we 
see from Eq. (7.1) that it is only the linear combina¬ 


tion ipo/fo + 4>il fi which dynamically relaxes against the 
strong CP angle 0. However, in the M —>■ oo decoupling 
limit, our model reduces to the standard single-axion 
model: the lighter mass eigenstate remains massless prior 
to the instanton-induced phase transition, gathers a mass 
A^aZ/o after this transition, and solves the strong CP 
problem. For this reason we shall continue to refer to 
the lighter mass eigenstate as our “QCD axion”, even in 
the presence of mixing. 

Given this, we now calculate the late-time energy den¬ 
sity associated with the lighter (QCD) axion as a func¬ 
tion of the lighter (QCD) axion mass x = Aq/oZmaj nor¬ 
malized to our decoupled benchmark limit with M = oo 
(for which a; = 1). (As a technical point, note that we 
describe this benchmark through the condition M = oo 
rather than the condition T = 1, since the former implies 
the latter regardless of fi/fo whereas the latter fails to 
imply the former when /iZ/o oo.) In all cases we take 
to = ^QCD and A(3 = Aqcd to have the specific fixed val¬ 
ues that correspond to the instanton-induced QCD phase 
transition, as plotted in Fig. 4. 

Our results are shown in Fig. 24. We see from Fig. 24 
that the late-time energy density associated with the 
QCD axion is suppressed relative to what would have 
been expected in the M —>■ oo decoupling limit. In¬ 
deed, this suppression is significant, occasionally stretch¬ 
ing over many orders of magnitude, and arises regard¬ 
less of the mass x of the QCD axion, its correspond¬ 
ing Peccei-Quinn scale /o, or the Peccei-Quinn scale /i 
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FIG. 24: The late-time energy density P)^^{Aqcd, M) associated with the QCD axion field <j)XQ, plotted as a function of the 
late-time axion mass x = Aq/o/ma S'lid normalized to what it would have been in the decoupling limit in which the second axion 
plays no role. These plots thus illustrate the effects of the non-zero mixing between our two axion fields in the presence of a 
non-zero QCD phase transition width Aqcd. Note that although each normalization factor p_;^|^(AQCD, M —>■ oo) is in principle 
different for each different /i//o curve, physics in the decoupling limit M —>■ oo is insensitive to the value of fi/fo- Thus all 
curves in each panel share the same normalization factor. 



XXX 

FIG. 25: Same as Fig. 24 except that the late-time energy densities (Aqcd, M) are now normalized to the values that they 
would have had if the QCD phase transition had been treated as instantaneous. These plots thus illustrate the effects that 
emerge due to the actual non-zero width of the QCD phase transition — effects which survive even in the M —>■ oo decoupling 
limit (for which a: —>■ 1). Note that unlike the plots in Fig. 24, here the normalization factors Pxg{0,M) are different for each 
/i//o curve and in fact even vary with x along each curve. 


associated with the second axion. Indeed, none of this 
would have occurred were it not for the non-zero width of 
the instanton-induced QCD phase transition. This then 
demonstrates the phenomenological relevance of this non¬ 
zero width, and the need to incorporate these kinds of 
effects in studies of multi-axion theories. 

Note that this suppression of the late-time energy den¬ 
sity due to the presence of the second axion field persists 
even in the limit for which /i 3> /q. However, this some¬ 
what counter-intuitive result should not come as a sur¬ 
prise. As we know, the limit fi/fo —f oo corresponds to 
one in which ^ f 0 — i.e., one in which the late-time 
mixing vanishes. However, this does not imply that the 
mixing is small throughout the time evolution — in fact, 
both 0 and its time derivative can be quite sizable dur¬ 


ing the mass-generating phase transition. Indeed, this 
is particularly true for the relatively small values of x 
which correspond to 77 < 0. Thus, it is during the mass¬ 
generating phase transition that the bulk of the energy 
density can be transferred to the second axion — even 
when fi fol This is therefore an effect that arises 
as a direct consequence of the non-zero width of the 
instanton-induced QCD phase transition and its associ¬ 
ated timescale. 

Interesting as these results are, they only address the 
issue of how our late-time energy density compares with 
what would arise in the decoupled case — i.e., the case 
with no mixing (^ = 0). In this paper, however, we have 
also repeatedly studied another comparison, namely to 
the case of an instantaneous phase transition — i.e., 
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the case in which we imagine our phase transition to 
have Aqcd = 0. The corresponding results are shown 
in Fig. 25. 

Once again, it is not difficult to understand the features 
in these plots. First, let us discuss the endpoints at x = 1. 
Of course, these endpoints correspond to the decoupling 
limit (because a; = 1 implies ^ = 0, as we have seen in 
Fig. 23); thus the physics at these endpoints includes the 
effects from only the lighter (QCD) axion itself. As we see 
from Fig. 25, even in this limit the resulting late-time en¬ 
ergy density is significantly suppressed by factors ranging 
from 5 to 500 as /o ranges from 10^^ GeV to 10^^ GeV. 
Thus the effects of properly accounting for the non-zero 
width Aqcd of our instanton-induced phase transition 
are extremely important, even for a single QCD axion! 

Moving away from the decoupling limit towards 
smaller values of x, we now begin to incorporate the ad¬ 
ditional effects of the mixing with the second axion field. 
As we see from Fig. 25, these effects enhance the result¬ 
ing late-time energy density above what it would have 
been in the decoupling limit. In so doing, they can even 
possibly overcome the above suppression completely — 
a prospect that depends on the value of fi/fo — and 
thereby potentially increase the late-time energy den¬ 
sity above what it would have been for an instantaneous 
phase transition! This enhancement ultimately reflects 
the same physics that we have already seen in Figs. 10 
and 11 for small mixing saturations in the to ~ re¬ 
gion. Indeed, we see from Fig. 23 that larger values of 
/i//o correspond to smaller mixing saturations 

Finally, we can also understand why these curves all 
begin at 1 for a: = 0. At this endpoint of these plots, 
our axion fields remain overdamped during the mass¬ 
generating phase transition. They thus remain insen¬ 
sitive to the transition width Aqcd- 

We can also consider the fraction of the total energy 
density of our two-axion system which is associated with 
the lighter axion field. These results are shown in Fig. 26. 
In general, we see that as fi/fo increases, an increasing 
fraction of the total energy density is associated with 
the heavier axion field. However, as a: —>■ 1, we see that 
virtually all of the total energy density can be associated 
with the lighter axion field. 

As a final comment, we remark that it is the combina¬ 
tion (j)o/fo + (f'llfi which not only solves the strong-CP 
problem but which also couples to gluons. Depending 
on the particular model under study, this combination 
may also couple to photons and other Standard-Model 
(SM) particles. We can therefore decompose this com¬ 
bination into our late-time mass eigenstates in order to 
determine the effective mass scales />„ and f\^ that gov¬ 
ern the relevant SM couplings of our mass eigenstates at 
any moment in time: 

^ ^ ^ ^ (7 7 ) 

fo fl ~ ho /Ai ■ ^ 

The late-time values of these effective coupling scales are 
plotted in Fig. 27. 
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FIG. 26: The late-time fraction Px^lp of the total late-time 
energy density which is associated with lighter axion field, 
plotted as a function of x. We see that this fraction generally 
decreases as /i//o increases, and also increases as a function 
of X. 


Several features of these plots are worthy of note. For 
example, we observe that in all cases the f^. coupling 
scales associated with our mass eigenstates are never sig¬ 
nificantly smaller than the fi coupling scales associated 
with our original unmixed states. Thus neither of our 
mass eigenstates will decay too rapidly to SM states, a 
situation which might potentially have resulted in phe¬ 
nomenological difficulty. Likewise, combining the results 
in Fig. 27 with those of Figs. 23 and 26, we see that 
there are a variety of different phenomenological situ¬ 
ations which can emerge, depending on the underlying 
parameters in our model. 

For example, let us consider the case with 
Mfo/^A ~ 0-4 and /i « 5/o. We immediately learn 
from Fig. 27 that our QGD axion coupling scale is 
/ao ~ 10^^ GeV, which is a full order of magnitude higher 
than the value that it would have had in the decou¬ 
pling limit. This implies that the QCD axion has an 
enhanced invisibility as assessed through its interactions 
with SM states. However, we see from Fig. 23 that these 
parameter combinations correspond to T « 0.4, where¬ 
upon we see from Fig. 26 that this in turn corresponds 
to an energy-density fraction Px^/p ~ 80%. This is thus 
a situation in which our QCD axion is invisible but nev¬ 
ertheless carries the bulk of the total late-time energy 
density. 

A somewhat different situation emerges if we keep 
fl ~ 5/o but now take M/o/^^ « 0.1. In this case, we 
learn from Fig. 27 that our QCD axion coupling scale 
is even higher — fxg « 10^^ GeV — which suppresses its 
couplings to SM states by three orders of magnitude com¬ 
pared with the decoupling limit. However, we see from 
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FIG. 27: The effective coupling scales fxg and f\^ associated with our axion mass eigenstates <j)Xo 4>Xi, plotted as functions 
of M/o//iA- We see that these scales are never significantly smaller than the fi coupling scales associated with our original 
unmixed states, and that either one or the other typically takes a value which significantly exceeds this when /i//o S> 1. 


Fig. 23 that these parameter combinations now corre¬ 
spond to 3; « 0.1, whereupon we see from Fig. 26 that 
this in turn corresponds to an energy-density fraction 
Pxg/p ~ 3%. This is thus a situation in which our QCD 
axion is invisible and carries almost none of the total 
energy density — all because of its mixing with the sec¬ 
ond axion in the presence of a non-zero phase transition 
width! 

The opposite situation, of course, emerges for large M. 
For example, with M/q/^a ~ 10 physics 

of this model is roughly insensitive to the value of /i//o; 
we find T « 1 and Px^/p ~ 100%. This, of course, sig¬ 
nifies nothing but the approach to our decoupling limit. 
However, as long as we are not precisely at the decou¬ 
pling limit, we still find from Fig. 24 that the energy 
density associated with this lighter axion field remains 
considerably suppressed compared with what we would 
have found at the full decoupling limit. Thus in this sense 
our QCD axion has regular couplings to SM states but 
nevertheless carries relatively little total energy density. 

Needless to say, our goal in this section has not been 
to propose a complete phenomenological model of ax¬ 
ion physics. Rather, it has merely been to illustrate the 
rich implications that can emerge for a QCD axion in 
the presence of mixing with a second axion when the 
non-zero width of the QCD phase transition is properly 
taken into account. As we have seen, the total energy 
density of such a two-axion system can be significantly 
suppressed relative to what might have been expected in 
the decoupling limit, and the individual energy densities 
associated with each axion can be severely distorted. In¬ 
deed, this lesson has been one of the primary themes of 
this paper. 


VIII. BEYOND TWO FIELDS 

In Sect. II, a toy model was constructed in order to 
study the effects that arise in multi-scalar models due to 
the simultaneous presence of two key ingredients: mass¬ 
generating phase transitions with finite widths, and non¬ 
zero mixing between the different components. For the 
sake of simplicity, our toy model consisted of only two 
scalar fields, as this is the minimum number that could 
be chosen for our purposes. However there is no reason 
that we are limited to two fields, and indeed the number 
of fields may be extended arbitrarily. Therefore, in order 
to gain a quick sense of what possibilities might emerge 
with an increased number of fields, we shall now briefly 
consider several aspects of the case with three fields. A 
more general study of the V-field case will be presented 
in Ref. [38]. 

In analogy with Eq. (2.8), we shall begin by assuming a 
rescaled (dimensionless) mass matrix Al^(t) which takes 


the form 

mloir) 


Wo2('r) 

M^{t) = 

m§i(r) 

1 -f to?i(t) 

771^2 (t) 



m?2('r) 

4-f TO^2( 


Here t = Mt (where M is the overall mass parameter 
which has been scaled out, as in Sect. H), and we assume 
that all components have a common time-dependence as¬ 
sociated with a mass-generating phase transition whose 
properties are exactly those of our two-component model. 
Note that the diagonal components of this matrix im¬ 
ply that our three fields have masses 0, M, and 2M at 
early times prior to the phase transition. These values 
are chosen for simplicity, and are also motivated by the 
Kaluza-Klein masses that might result in a theory with 
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FIG. 28: The region within the three-dimensional mixing 
space parametrized by {fnQi,rnQ2,7ni2} for which the mass 
matrix in Eq. (8.1) is positive-semidefinite. For this plot 
we have taken thoo = 5, rhii = 5, and m 22 = T 



FIG. 29: The two smaller mass eigenvalues Aq and Xf in onr 
three-field model, plotted as a function of time for different 
values of mo 2 = hli 2 - As in Fig. 28, we have taken ttiqo = 5, 
mil = 5, and 77122 = T We see that each field can now sepa¬ 
rately have its own pulse with its own effective frequency, re¬ 
sulting in the possibility of two parametric resonances which 
can either reinforce or interfere with each other. 


a flat extra dimension (in which case we would identify 
M as the inverse length of this dimension). 

Inherent in this mass matrix are now three independent 
mixing angles which correspond to rrigi, mg 2 , and mf 2 - 
However, just as in the two-field case in Sect. II, there 
are ultimately constraints on these mixing parameters 
which ensure that our mass matrix remains positive- 
semidefinite. Such constraints are illustrated in Fig. 28 
for the case with m^g = 5, fhii = 5, and rh 22 = 1 — a 
fixed parameter choice we shall adopt for all of the plots 
in this section. 

In our two-component toy model, we have seen that the 
lighter field can occasionally develop a “pulse” which is 
ultimately responsible for the parametric-resonance and 
re-overdamping phenomena in Sects. V and VI respec¬ 
tively. However, in a t/iree-component model, it turns 
out that the two lightest fields can each develop a pulse. 
This is illustrated in Fig. 29. Indeed these two pulses 
need not have the same effective frequencies or corre¬ 
spond to the same values of Aq- As we shall see, this can 
then give rise to multiple resonances which can reinforce 
or interfere with each other and thereby produce new ef¬ 
fects which transcend those realizable with only a single 
pulse. We also note that the pulse of our second mass 
eigenvalue can leave this eigenvalue with a smaller 
magnitude at late times than at early times — a phe¬ 
nomenon which was not possible for Ag in the two-field 
case. This has the effect of allowing a further elongation 
of the time interval during which this field can remain in 
a re-overdamped or parametrically resonant state. 


To illustrate the interplay between the different res¬ 
onances, we can consider the late-time energy densities 
of this system. Since we now have three individual con¬ 
tributions to the total energy density, we shall highlight 
their separate effects by considering the partial sums 

i 

Pi = ( 8 . 2 ) 

i=0 

In Fig. 30 we have plotted these partial sums, with 
each normalized to the value that it would have had if 
the mass-generating phase transition had been instanta¬ 
neous. This demonstrates the combined effects of the 
mixing and the finite phase-transition width. The left 
panel of Fig. 30 illustrates the behavior of Pq = Pao 
alone: this behavior is similar to what we have already 
seen in our two-field toy model, for which the pulse 
in the evolution of Ag produces resonant peaks in the 
corresponding energy-density component at late times. 
These peaks grow stronger for larger TOq 2 , with the en¬ 
ergy density in this field dissipating inversely with Ag 
once Aq > 27r/Ag. However, the center panel shows the 
combined contribution to the total energy density from 
the two lighter fields and thus exhibits something new: 
the contributions of the two separate parametric reso¬ 
nances combine to reinforce each other and even produce 
an effective enhancement plateau which stretches over an 
extended interval in Aq. Indeed, all along this plateau, 
the energy densities of our system can be enhanced by 
as much as nearly an order of magnitude! Finally, the 
right panel of Fig. 30 illustrates the behavior of the to- 
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FIG. 30: The partial-sum energy densities Pt{Ao), i = 0,1,2, plotted as functions of the phase-transition width Aq and 
normalized to the values that they would have had for an instantaneous phase transition. As in Figs. 28 and 29, we have 

again taken moo = 5, mn = 5, and m 22 = 1. We observe many features familiar from our two-component toy model, including 
parametric resonances which enhance the corresponding energy densities as well as suppressions by many orders of magnitude 
that develop for sufficiently large Aq. However, there are also new features which arise due to the fact that we now have more 
than two fields in our system. The center panel, in particular, illustrates the manner in which the separate resonances from 
the two lighter fields reinforce each other and produce a plateau along which the partial-sum energy density Pi is enhanced by 
nearly an order of magnitude over a significant A( 3 -interval. 


tal energy density p. Note that in this case, the heaviest 
field effectively dominates the energy density for < 5, 
thereby washing out those parametric resonances that 
appear in that range, while leaving those that appear for 
Ag>5. 

The purpose of this section has merely been to pro¬ 
vide several short examples of new phenomena that can 
arise when more than two scalar fields are involved. We 
shall defer a more detailed study of these general cases 
to Ref. [38]. 


IX. CONCLUSIONS 

In this paper we have examined some of the novel ef¬ 
fects that emerge when multiple scalars mix with each 
other and experience a time-dependent mass-generating 
phase transition within the context of a time-dependent 
background cosmology. As we have explicitly demon¬ 
strated within a relatively simple two-field toy model, 
these effects include both enhancements and suppressions 
of the total late-time energy density of the system as well 
as the late-time energy densities of the individual scalar 
components. These enhancements and suppressions need 
not be mere 0{l) effects; indeed, we have seen that the 
resulting late-time energy densities can often be altered 
by many orders of magnitude. These alterations are even 
sufficient to produce wholesale shifts in the identity of 
the field carrying the bulk of the late-time cosmological 
abundance. 

We also found that our late-time energy densities can 
exhibit large parametric resonances which make them ex¬ 
tremely sensitive to relatively small changes in the degree 
of mixing or in the width of the phase transition. Thus, 


any effects which might change — even slightly — the 
rate at which the phase transition unfolds can have a 
significant impact on the late-time abundances of our 
fields. We also demonstrated that our system can en¬ 
ter a so-called “re-overdamped” phase in which the field 
values and energy densities behave in a manner quite 
different from what would normally be associated with 
pure dark matter or vacuum energy. Indeed, we saw that 
the system can remain in this re-overdamped phase for a 
considerable length of time, suggesting that this phase be 
considered alongside “vacuum energy” and “matter” as 
an equally valid new behavioral phase that a scalar field 
can experience during a significant portion of cosmolog¬ 
ical evolution. Finally, we illustrated the importance of 
some of these effects within the context of a simple model 
in which a second axion is incorporated into the standard 
QCD axion theory. 

As indicated in the Introduction, our results can be 
taken as an explicit demonstration of how critical it 
can be to properly incorporate the non-zero but finite 
widths of mass-generating phase transitions when dis¬ 
cussing late-time energy densities in cosmological settings 
involving multiple mixed scalars. Indeed, treating these 
widths as being either zero or infinite (the latter corre¬ 
sponding to a so-called “adiabatic” approximation in the 
axion case) can lead to results for late-time quantities 
such as cosmological abundances which are significantly 
distorted from their true values, often by many orders of 
magnitude! Given the sensitivity with which the differ¬ 
ent slices of the so-called “cosmic pie” have already been 
measured, such effects will inevitably play critical roles in 
determining whether a particular cosmological scenario is 
viable or ruled out. 

Needless to say, the features we have discussed in this 
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paper give rise to many new possibilities for the phe¬ 
nomenology of the early universe and cosmological evo¬ 
lution. As such, new approaches to model-building can 
readily be imagined. For example, many extensions to 
the Standard Model involve scalar fields which are mass¬ 
less (or effectively massless) at high scales, but for which 
masses are generated dynamically at lower scales. Ex¬ 
amples include string moduli, string axions, and other 
axion-like particles. However, these particles can be cos- 
mologically problematic if these dynamically-generated 
masses are sufficiently large because they can overdose 
the universe once they begin oscillating and act as mas¬ 
sive matter. Indeed, this is the crux of the cosmological 
moduli problem [41-43]. Moreover, particles of this sort 
are often produced through mechanisms such as vacuum 
misalignment which are completely decoupled from the 
dynamics of inflation. In such cases, spatial fluctuations 
in the energy density of these particles therefore give rise 
to isocurvature perturbations, the observational limits on 
which are quite severe [44]. 

Fortunately, mechanisms such as those we have dis¬ 
cussed in this paper can suppress the resulting late-time 
energy densities by many orders of magnitude. This en¬ 
hanced energy dissipation can therefore provide a novel 
way of alleviating such constraints on models of new 
physics, and may represent an alternative solution to the 
moduli and isocurvature problems. They may also allow 
a general weakening of the cosmological bounds on the 
energy scales associated with the new physics (such as 
the Peccei-Quinn scale associated with axions). 

As we have discussed, the features explored in this pa¬ 
per can lead not only to a reduction in the total late¬ 
time energy density of our system, but also to changes 
in the late-time partitioning of that total energy density 
amongst the fields in our theory. In new physics sce¬ 
narios which involve multiple scalar fields, this too can 
have observable phenomenological consequences. For ex¬ 
ample, certain classes of string models are generically 
expected to give rise to large numbers of light axion-like 
particles [30-33] , the masses of which must be generated 
dynamically via non-perturbative effects. In situations 
in which the same dynamics contributes to the genera¬ 
tion of mass terms for multiple such particles, the relative 
cosmological abundances of those particles can be signif¬ 
icantly altered by the effects we have studied here. 

There are many ways in which the analysis in this pa¬ 
per might be extended and generalized. For example, 
in this paper we have focused primarily on scalar fields 
which have a negligible back-reaction on the background 
cosmology, at least during the epoch of dynamical mass 
generation. In other words, we have assumed that our 
scalars do not contribute significantly to the total energy 
density of the universe until well after their masses have 
settled into their late-time values. However, in principle, 
we could relax this assumption and broaden this study 
to include scalars which do have a significant impact on 
the background cosmology. This would then establish a 
non-linear back-reaction which could possibly enrich the 


dynamics of this system even further, with the evolution 
of the Hubble parameter H{t) itself depending on the 
behaviors of the scalar energy densities as functions of 
time. 

This possible extension of our work could be partic¬ 
ularly important if the scalars involved are those which 
are directly responsible for triggering periods of rapid 
cosmological expansion, including cosmic inflation. In¬ 
deed, a wide variety of inflationary models involving mul¬ 
tiple scalar fields exist in the literature; these include 
hybrid inflation [45, 46], assisted inflation [47-50], N- 
flation [51], and multi-field stochastic inflation [52]. Like¬ 
wise, several constructions have been proposed through 
which a potential for the fields responsible for cosmic in¬ 
flation is generated dynamically [53-56]. Such construc¬ 
tions are typically motivated by the possibility of real¬ 
izing a viable inflationary potential without the intro¬ 
duction of extremely small couplings or arbitrary mass 
scales. It would therefore be extremely interesting to ex¬ 
plore the cosmological consequences of a time-dependent, 
dynamically-generated inflaton potential in the context 
of multi-field inflation. 

In this connection, it is perhaps worth emphasiz¬ 
ing that the parametric resonance we have discussed in 
Sect. V is entirely unrelated to the parametric resonance 
which gives rise to an epoch of explosive production 
of light scalars in inflationary models with a preheat¬ 
ing phase [57, 58]. Indeed, the effect we have discussed 
in this paper involves a resonance between the oscilla¬ 
tion frequency of the light scalar, as given by its mass 
Aq, and the effective modulation frequency uJeS of that 
mass. As such, this parametric resonance is ultimately 
driven by the dynamics of the mass-generation mecha¬ 
nism and leads to an enhancement in the amplitude for 
coherent oscillations of the zero-momentum mode of (j)\g. 
By contrast, the effect discussed in preheating scenar¬ 
ios involves a resonance between the frequency of the 
coherently-oscillating inflaton field and the frequency as¬ 
sociated with a given momentum mode of the lighter 
scalar field into which it decays. Such a resonance is 
ultimately driven by inflaton dynamics and results in 
particle-like excitations of this lighter scalar. 

Other extensions and generalizations of our work are 
also possible. For example, our original dynamical equa¬ 
tions in Eq. (2.1) assume that there are no significant 
spatial variations in our field values, but in general it is 
possible to allow our field values to vary not only in time 
but also in space. This could then give rise to distinctly 
different physics in different spatial regions, with domain 
walls potentially separating different regions experienc¬ 
ing different phases and different energy densities. Such 
spatial variations and domain walls would inevitably in¬ 
troduce additional length scales into the problem, and 
thereby likely give rise to new complex dynamical be¬ 
haviors for our system as a whole. 

In a similar vein, we could likewise imagine that our 
fundamental dynamical equations in Eq. (2.1) are modi¬ 
fied by the inclusion of source terms reflecting possible in- 
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teractions with more complex cosmological environments 
{e.g.^ with other fields). Such source terms could then 
give rise to additional features which effectively modify 
our mass eigenvalues and eigenstates, trigger enhance¬ 
ments or suppressions of late-time energy densities, and 
produce additional kinds of resonances. Indeed, such 
effects would be analogous to the well-known “matter 
effects” {i.e., MSW effects) that arise when neutrinos 
propagate through matter rather than through a pure 
vacuum. 

Another possible generalization of our toy model is to 
consider more than two scalar fields. Indeed, as briefly 
discussed in Sect. VIII, there is no limit to the number 
of scalar fields which may be considered. These models 
can therefore grow in both complexity and sophistica¬ 
tion compared with the two-component model we have 
studied here. 

One natural possibility along these lines is to con¬ 
sider the case in which the different fields (j)i are noth¬ 
ing but the Kaluza-Klein (KK) modes of a single higher¬ 
dimensional field <&. This is indeed a well-motivated 
possibility, as string theory naturally gives rise to many 
scalars (axions, geometric/gravitational moduli, gauge- 
neutral singlets, and so forth) that populate the “bulk” 
volume transverse to the brane on which gauge inter¬ 
actions reside. As such, the different KK modes of 
such bulk fields are natural dark-matter candidates, even 
though they might be unstable — an observation which 
has motivated the Dynamical Dark Matter framework 
of Refs. [28, 29, 36]. In this case, the masses that we 
have assumed to exist prior to the phase transition would 
be nothing but the corresponding KK masses, while 
the extra time-dependent contributions mijh^r) are pre¬ 
sumably generated by a phase transition which breaks 
the higher-dimensional symmetries that would otherwise 
have aligned our mass eigenstates with KK eigenstates. 
The resulting mixing structure is therefore determined 
according to the intrinsically higher-dimensional geomet¬ 
ric symmetries of the theory and the manner in which 
these symmetries might be broken by the phase transi¬ 
tion. (For example, in the five-dimensional axion scenario 
of Refs. [29, 35, 36] the instanton-induced phase transi¬ 
tion takes place purely on the brane and thereby breaks 
a translation symmetry that otherwise existed along the 
extra dimension.) In general, the cosmological properties 
of such KK systems can be extremely rich, and will be 
discussed in more detail in Ref. [38]. 

Given these observations, it might at first glance seem 
that the two-component toy model we have studied in 
this paper — while illustrative for pedagogical purposes 
— might not have been worthy of the detailed attention 
we have afforded it. However, this model can always be 
viewed as the low-energy limit of a more complete model 
(such as the model of Kaluza-Klein cosmology discussed 
above) in which only the contributions from the lightest 
two modes are considered. Indeed, in many cosmolog¬ 
ical settings, the contributions of the lightest fields are 
likely to dominate the resulting phenomenology and/or 


carry the largest abundances. Thus, in this sense, our 
two-component toy model can be viewed as the common 
root underlying a variety of much more general models 
— a “kernel” which gives rise to rather universal features 
and behaviors that will continue to appear (perhaps with 
further embellishments) in more complete or realistic set¬ 
tings. As such, then, this toy model can be viewed as rep¬ 
resenting a universal low-energy limit of a wide class of 
models of this type whose leading-order phenomenologies 
will be exactly those we have investigated in this paper. 
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Appendix A: Narrow-width behavior of the 
late-time energy density: Analytical results 

As we have seen, the presence of a non-zero width Aq 
for the mass-generating phase transition prevents us from 
obtaining analytical results for most of our quantities of 
interest. However, in the case with Ac = 0 (correspond¬ 
ing to an instantaneous phase transition), the resulting 
dynamics can be treated analytically. In this Appendix 
we provide analytical results for the late-time energy den¬ 
sity p. 

For Ac = 0, the behavior of each of our two fields is 
relatively simple. Prior to the phase transition at tq, 
only the massless field has a non-zero value Aq, as in 
Eq. (2.16); this field is necessarily overdamped, however, 
and thus remains effectively constant until the phase 
transition at tq. At that point, the mass matrix receives 
the extra contribution which mixes the two fields with 
an angle 0; the initial massless field value is thus redis¬ 
tributed at T = TG such that 

'/ao (tg) = ^0 cos 0 , (/>Ai (tg) = Ao sin 0 . (Al) 

Each field then evolves forward in time independently. 

At first glance, motivated by the behavior illustrated 
in Eig. 1, we might try to estimate the total late-time 
energy density p using a virial approximation in which we 
assume that the energy density associated with each field 
(p\. remains constant until = Kj2Xi [where k = 2 (or 
3/2) for a matter- (radiation-) dominated universe], after 
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which it scales as ~ r Defining t\ = niax{rG, k/2X}, 
we would then write 

^ ^ E (^) = ^ E iMrc)? . 

' (A2) 

Unfortunately, this approximation breaks down for any 
fields 4>i for which the threshold of critical damping 
= nj^Xi is close to the phase transition time tq. 

For this reason, we shall instead evaluate our field val¬ 
ues and corresponding energy densities exactly. In the 
limit of an instantaneous phase transition occurring at 
T( 5 , the exact solution for each field is given by 

<^A(r) = I [(j^xiTa)] >^tg X 

[Jk+ iXTG)Y^_ (At) - (Atg) Jk_ (At)] 

(A3) 

where k± = 1)/2, where A refer to the mass eigen¬ 

values for T > TG, and where J,y{z) and IU(z) are Bessel 
functions of the first and second kind, respectively. From 
this result we may directly calculate the corresponding 
energy densities for all t >tg, obtaining 

P\{t) = y [<(>a(tg)]^ A^(tgt) Jxir) (A4) 

where 

Jx = [Jn+ iXTG)Yi^_ (At) - {Xtg)Jk_ (At)] ^ 

+ [•/«+ (Atg)Fk4. (At) (Atg)</k^ (At)] 

(A5) 

Given these exact results, we can now extract the late¬ 
time behavior of p\ through the asymptotic expansions 
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which hold for z ^ 1. With this approximation we find 
that 


Jx 


— [J^,(Atg)+U4(Atg) 


(A7) 


whereupon use of Eq. (Al) yields the total late-time en¬ 
ergy density 


a =\aI irar) (i) x 

{Xt,Jx„cos'^e + X\Jx,syn^e) . (A8) 

The result in Eq. (A8) is an exact expression for the 
late-time energy density p in the Ag —>■ 0 limit. Given 
this result, several things are immediately clear. First, 


we observe that in the case of matter-dominated universe 
(corresponding to k = 2 ), the total late-time energy den¬ 
sity simplifies to 

2 + (Aq + Ai)tg -I- (Ag — X\)tg cos(20) 

= £1(1 + <Jg) , (A9) 

where in passing to the second line we have used the re¬ 
sults in Eq. (2.32). As a result, all 0-dependence drops 
out of the final result for the total late-time energy den¬ 
sity. In other words, for a matter-dominated universe, the 
corresponding late-time energy density is independent of 
mixing] 

Second, we see that the virial approximation in 
Eq. (A2) corresponds to the limit in which tgA^ ^ 1 for 
both fields. Indeed, in this limit we see that t\ = tg 
for each field (since both fields are already significantly 
underdamped when the phase transition occurs), and we 
can likewise apply the asymptotic expansions in Eq. (A6) 
to the remaining Bessel functions in Eq. (A7), thereby 
reproducing the result in Eq. (A2). Further use of the 
results in Eq. (2.32) then gives the total late-time energy 
density 



(AlO) 

which is again independent of the mixing — now regard¬ 
less of the value of k\ However, comparing this result 
with the exact result in Eq. (A9) for the case of a matter- 
dominated universe illustrates the imperfect nature of the 
virial approximation. 



Appendix B: Late-time energy densities: An 
alternative approach 

In the main body of this paper, we have consistently 
viewed the matrix elements as our independent vari¬ 
ables. These have been reparametrized in other forms, 
giving rise to variables such as Am^, rj, 9, and 

but in all cases these independent variables have served 
as inputs. We then calculated the values of dependent 
quantities such as the eigenvalues A^, energy densities p, 
and so forth. Thus, when we have studied how a certain 
dependent variable such as a late-time energy density p 
depends on the mixing, we have implicitly held the diag¬ 
onal matrix elements fugg and rhii (or equivalently 
and p) fixed, and varied the off-diagonal matrix element 
fhQi (or equivalently the mixing angle 9 or mixing satu¬ 
ration ^). 

However, another possibility is to study the mixing- 
dependence of quantities such as p when we hold the 
late-time eigenvalues Xi fixed instead. This approach 
can be useful in situations where the masses of our states 
are known and it is only the mixing between such states 
that we wish to vary. 
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In order to implement this approach, we must first in¬ 
vert our usual algebraic relations and express {itiqq, mh} 
in terms of {ApjAf}. However, we must also recognize 
that our previous assumption of holding {mpg, rriii} fixed 
is also implicitly buried within the definition of another 
variable: the mixing saturation Although uiqi or 9 
parametrize the mixing in absolute terms and are thus 
independent of such assumptions, the mixing saturation 
^ defined in Eq. (2.28) describes the mixing as a frac¬ 
tion of the maximum allowed mixing, and this maximum 
allowed mixing has always been determined relative to 
a fixed {inQQ,m1i}. What we now need, by contrast, is 
a new variable which describes the mixing as a frac¬ 
tion of the maximum mixing allowed for a given fixed 
{Ag, Aj}. Indeed, this is an entirely different quantity. 

Thus, we seek to express {mQQ,mli,^} in terms of 
{Ag,Af,^A}- Equivalently, defining 


and 


= -^0 + "^1 

AA^ ^ Af-Ag 

"^0 + "^1 


(Bl) 


(B2) 


in analogy with and 77 , we would like to express 

in terms of Surprisingly, 

these relations ultimately take a relatively simple form: 


^sum ^sum 


2 2 
V = Vx 


^-ex 

C = ■ 

These relations in turn imply that 


m, 


01 




(B3) 


(B4) 


Likewise, following Eq. (2.30), we can also determine the 
absolute mixing angle 6 via 

tan2(20) = ■ (B5) 

Note that 0 < ? 7 a < 1, whereas —1 < rj < 1 . Also note 
that the eigenvalues Ag are independent of the sign of 77 
(as evident from Eig. 6 ). It is for this reason that our 
inverse relations in Eq. (B3) determine 77 [and likewise 
tan(20)] only up to an overall sign. Thus, strictly speak¬ 
ing, our variable map really trades { 77 , C} for {77A,CA,e} 
where e = sign(? 7 ) = ±1 is a discrete Z 2 phase. Having 
e > 0 corresponds to 77 > 0 and tan(20) > 0. However, 
we shall avoid writing e in what follows, understanding 
the proper phase choice to be implicit. 

In principle, given this variable map, we could now 
work our way through this entire paper again using 
the variables {Ag^^, 77 A, Ca} rather than the variables 
Indeed, every curve in every plot which 


has been labeled in terms of certain numerical values 
for the first set of parameters could now equivalently be 
relabeled in terms of different numerical values for the 
second set of parameters. In cases where these curves 
describe the behavior of a certain quantity such as an 
energy density as a function of time, or the behavior of 
a quantity such as a late-time energy density as a func¬ 
tion of a phase-transition width Ac, these curves would 
themselves remain intact and continue to take the same 
shapes as before because our time and Aq parameters 
are independent of our variable change. All that would 
change are the labels associated with each such curve. 
Indeed, every feature that appears in our plots for cer¬ 
tain fixed values of {msumiV} (such as the oscillations 

due to parametric resonances or the re-overdamping of 

_2 

certain fields) would continue to exist when {Ag^^, 77 ^;,,} 
are held fixed instead. As a result, nothing we have 
said in this paper concerning the parametric-resonance 
or the re-overdamping phenomena would be affected by 
this change in variables. 

The only changes that could potentially occur due to 
this variable exchange concern the effects of mixing as 
discussed in Sects. HI and IV. Of course, each individ¬ 
ual curve within the plots in Figs. 10, 11, 12, 16, and 18 
would remain exactly as before and would merely be rela¬ 
beled using {Agmjj, 77 a, ^a} rather than 77 , ^}. How¬ 

ever, what would change is the resulting grouping of these 
curves into distinct panels in which all curves share com¬ 
mon values of, say, ijx rather than 77 . Indeed, although the 
space of all possible curves is not changed, these curves 
would experience a reshuffling into differently-grouped 
subsets. As a result, the only physics questions whose 
answers might change in this process are those which 
rely on comparisons between the different curves within 
a given grouping. Clearly, questions pertaining to the ef¬ 
fects of mixing (ie., the effects of changing the values of 
^ or ^A while holding 77 or px fixed) are in this category. 

Rather than provide an exhaustive study of such re¬ 
groupings and their outcomes, we shall here content our¬ 
selves with providing a concrete example of the effects 
that can emerge by focusing on a single pair of curves 
which we shall call A and B. These curves can be taken 
to show, for example, a late-time energy density p as 
a function of Aq, but they will differ in their absolute 
mixing by an angle Ad when are held fixed. 

We will then calculate a third curve B' which also differs 

from the curve A by the same absolute angle A9 when 
_2 

{Agumj^A} fixed instead. In this way, comparing 

curves B and B' will then give us an idea of the difference 
that can come from the choice of which variables to hold 
fixed. 

For concreteness, we select curves A and B to be 
the ^ = 0.9 and ^ = 0.5 curves within the right panel 
of Fig. 11. Within this panel, Afn^ = 4 and 77 = 0.75 

for both curves. Use of Eq. (B3) then provides us 

_2 — 

with the corresponding {Ag^JJJ, values for each of 

these curves: we find {4.0,0.96, 0.94} for curve A and 
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FIG. 31: The late-time energy density p, plotted as a function of Aq in a matter-dominated nniverse (left panel) or radiation- 
dominated nniverse (right panel) and normalized to its value for an instantaneons phase transition with zero mixing. In each 
case, the curves A, B, and B' correspond to the parameter choices in Table I: curves A and B differ by a fixed A9 = 6b — 9a 
while holding fj fixed, while curves A and B' differ by the same A9 but instead hold fixed. Thus, the difference between 
curves B and B' in each case illustrates the impact of holding p fixed versus holding p^ fixed when performing a shift of the 
mixing angle A9 relative to the common curve A. 


{4.0,0.82,0.61} for curve B. Likewise, use of Eq. (B5) 
with either the jly,or {??>,,^>,} variables then tells us 
that 9a ~ 19.22° and 9b ~ 11.90°. [The fact that these 
values emerge in each case using either set of variables 
is a useful cross-check that the absolute mixing angle 
9 is indeed invariant across the mapping in Eq. (B3).] 
We now wish to define a third curve B' such that px is 
held fixed relative to curve A (be., p^^ ^ = p^x'^ « 0.96) 
while the mixing-angle shift A9ab' = 9b' ~ 0a matches 
A9ab = 9b — 9a- This latter condition implies that 
9b' =9b~ 11.90°, whereupon use of Eq. (B5) tells us 

that ^^x ^ ~ 0.84. Use of Eq. (B3) then tells us that 
~ {4, 0.88, 0.80} for curve B'. These param¬ 
eter values for curves A, B, and B' are summarized in 
Table I. 



P 


Vx 

Ca 

9 

A 

0.75 

0.90 

0.96 

0.94 

19.22° 

B 

0.75 

0.50 

0.82 

0.61 

11.90° 

B' 

0.88 

0.80 

0.96 

0.84 

11.90° 


TABLE I: Parameters defining the curves A, B, and B' shown 
in Fig. 31. All curves correspond to fn^uni = -^sum = 4 and 
e = -t. 

In Fig. 31 we plot Curves A, B, and B'. The left panel 
of Fig. 31 shows these curves for a matter-dominated 
universe, while the right panel shows these same curves 
for a radiation-dominated universe. Thus, while curves A 
and B in the left panel of Fig. 31 respectively correspond 
to the ^ = 0.9 and ^ = 0.5 curves in the right panel of 
Fig. 11, curves A and B in the right panel of Fig. 31 


respectively correspond to the ^ = 0.9 and ^ = 0.5 curves 
in the right panel of Fig. 10. 

As evident from Fig. 31, there are significant differ¬ 
ences between curves B and B'. Indeed, these curves 
are sufficiently distinct in the radiation-dominated case 
(right panel) to turn what would have been an overall 
enhancement of the normalized late-time energy density 
(curve B) into an overall suppression (curve B')\ We 
thus conclude that the question of which variables to 
hold fixed can indeed be a critical one for certain types 
of physics questions within certain regions of parame¬ 
ter space. This also highlights the need to be extremely 
careful about specifying which parameters are held fixed 
when making statements concerning the effects of vary¬ 
ing some parameters relative to others. In this paper, 
in order to study the effects of mixing, we have consis¬ 
tently held our diagonal mass-matrix elements mgQ and 
mil fixed while varying m^i (or equivalently varying 9 or 
^). This ultimately represents a choice which stems from 
our desire to perturb the potential generated by our un¬ 
derlying mass-generating phase transition as minimally 
as possible. However, this issue will be explored further 
in Ref. [38]. 
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